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ABSTRACT 


This  report  presents  a  user-friendly  Interactive  computer  program, 
written  In  Fortran  IV,  that  tabulates  and  evaluates  the  single  and  double 
sampling  plans  contained  In  MIL-STD-105D.  The  program  also  provides  the  user 
an  option  to  plot  curves  for  the  associated  operating  characteristics  when  the 
switching  rules  are  applied. 

The  sampling  scheme  is  described  in  terms  of  a  Markov  chain  to  obtain 
steady-state  probabilities  of  being  in  the  various  states  of  normal,  tightened 
and  reduced  inspection.  From  these  probabilities,  such  measures  of 
effectiveness  as  the  operating  characteristic  curve,  average  sample  number, 
average  outgoing  quality,  and  average  fraction  inspected  are  obtained^- 


ii 
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CHAPTER  1 
INTRODUCTION 

The  ABC-STD-105  represents  the  culmination  of  an  effort  made  by  a  working 
group  comprised  of  the  military  agencies  of  the  USA,  Great  Britain  and  Canada 
In  developing  a  common  standard  for  acceptance  sampling  by  attributes.  The 
complete  procedures  and  tables  were  approved  and  published  by  the  U.S. 
Department  of  Defense,  Washington,  D.C.,  as  Military  Standard  105D  (MIL-STD- 
105D)  on  April  29,  1963.  [181.  Henceforth  MIL-STD-105D  will  be  referred  as 
the  Standard. 

The  Standard  Is  applicable  to  a  continuing  series  of  lots  where  a  lot 
sometimes  Is  defined  as  all  of  the  material  received  at  particular  time.  In 
other  cases,  the  lot  might  be  all  of  the  product  manufactured  over  a  specific 
time  period  .  The  aim  in  defining  a  lot  is  to  obtain  a  homogeneous  grouping 
of  products.  It  is  desirable  to  have  large  lots  to  Improve  the  scheme's 
ability  to  discriminate  between  good  lots  and  bad  lots  for  comparable  total 
sampling  effort. 

The  usefulness  of  any  sampling  plan  can  be  measured  by  how  well  the  plan 
Identifies  poor  quality  product,  and  exerts  pressure  upon  the  producer  to  take 
action  to  Improve;  MIL-STD-105D  Is  intended  to  apply  such  pressure  on  the 
producer  of  a  continuing  supply  of  product  to  produce  at  a  level  of  quality 
equal  to  or  better  than  a  designated  Acceptable  Quality  Level  (AQL).  The 
Standard  defines  the  AQL  as 

.  .  .  the  maximum  percent  of  defective  (or  the  maximum 
number  of  defects  per  hundred  units)  that,  for  purposes  of 
sampling  inspection,  can  be  considered  satisfactory  as  a 
process  average.  [18:31 

To  exert  such  pressure,  a  set  of  rules  for  switching  from  normal  to 
tightened  Inspection  and  back,  based  upon  the  results  of  samples  taken  from 


recently  submitted  lots,  is  an  integral  and  necessary  part  of  the  Standard. 

In  addition,  a  reduction  to  reduced  sampling  following  favorable  experience  is 
also  Included.  It  is  the  incorporation  of  these  switching  rules  which  makes 
the  standard  a  sampling  system. 

The  procedure  for  switching  between  plans  is  Important  since  It  is 
designed  so  as  to  effect  optimum  protection,  both  to  its  producer  as  well  as 
to  the  consumer.  When  the  quality  of  the  produced  goods  falls,  corrective 
action  is  instituted  and,  if  the  quality  Improves,  there  is  a  reward  in  the 
shape  of  a  less  severe  inspection  plan. 

The  MIL-STD-105D,  which  has  become  a  standard  for  lot-by-lot  sampling 
inspection  by  attributes  for  all  Industry,  it  not  only  applicable  to  "nuts  and 
bolts"  but  also  the  sampling  inspection  of  parts,  components,  subassemblies, 
assemblies  and  end  items.  It  also  may  be  used  in  administrative  situations, 
such  as  in  audit  sampling  for  record  keeping  accuracy,  accounts  receivable, 
workload,  inventory  and  so  forth. 

The  importance  of  MIL-STD-105D  can  best  be  characterized  by  a  statement 
in  a  paper  edited  by  William  R.  Pabst,  Jr.,  which  is  as  follows: 

In  the  years  of  usage  of  1050,  no  serious  complaint  and 
no  pressure  for  revision  have  been  evidenced  by  those  using 
the  standard,  either  in  the  structure,  the  detail  or  the 
clarity.  .  .  .  MIL-STD-105D  Is  the  most  widely  used  of  all 
the  acceptance  sampling  schemes.  [15:99] 

Objective 

The  main  purpose  this  study  was  to  provide  a  tool  for  the  evaluation  of 
the  MIL-STD-1050  system  of  sampling  plans  by  the  design  and  implementation  of 
a  completely  interactive  FORTRAN  IV  language  computer  program. 

Three  main  objectives  were  emphasized.  First,  the  sample  size  code 
letter,  sample  size,  acceptance  and  rejection  number  data  as  specified  in 


Tables  I,  II-A,  8,  C,  III-A,  B,  C  of  the  Standard  for  single  and  double 
sampling  plans  are  made  available.  Thus  users  can  obtain  the  “raw"  data  that 
would  be  needed  for  implementation  and  evaluation  of  the  Standard.  Second, 
once  the  plan  parameters  are  obtained  for  normal,  tightened  and  reduced 
inspection  levels  (at  the  designated  AQL),  users  have  the  option  of  analyzing 
the  system  and  have  at  their  disposal  various  measures  of  effectiveness 
including  tabulations  of  probability  of  acceptance  ( OC  Curve)  values.  Average 
Sample  Size  (ASN),  Average  Outgoing  Quality  (A00)  values  and  the  Average 
Fraction  Inspected  (AFI)  for  any  given  fraction  defective  values  input  to  the 
program. 

Third,  the  computer  program  can  provide  plots  of  the  operating 
characteristic  (OC),  AOQ,  ASN  and  AFI  curves  for  the  plans  that  are  under 
scrutiny. 

Related  Work 

MIL-STD-105D  is  a  sampling  scheme  the  correct  implementation  of  which 
combines  several  individual  sampling  plans  into  a  procedure  designed  to  use 
economic,  psychological  and  operational  means  to  motivate  a  supplier  to 
sustain  quality  at  least  at  the  level  of  a  prescribed  AQL.  Btrt~th1s 
motivation  is  not  restricted  to  the  supplier  alone.  Since  the  revision  of  the 
Standard  by  the  ABC  Working  Group  and  its  publication  by  the  U.S.  government 
in  1963,  several  research  leaders  In  the  areas  of  quality  control  and 
statistics  took  up  the  task  of  analyzing  and  evaluating  the  complexities 
Involved  with  the  operation  of  this  sophisticated  standard. 

Use  of  any  specific  sampling  plan  without  utilizing  the  switching  rules 
Is  considered  nothing  less  than  the  misuse  of  the  Standard  itself. 
Nevertheless,  the  Standard  presents  OC  curves,  values  of  AOOL  and  limiting 
quality  (LQ)  for  specific  plans  only  and  not  for  the  system  as  a  whole. 


Limitations  at  the  time  of  publication  in  the  availability  of  computational 
procedures  to  determine  the  operating  characteristics  and  other  effectiveness 
measures  of  the  system  was  mostly  responsible  for  that  situation. 

In  1965,  H.  F.  Dodge  [3]  presented  a  method  for  determining  the  operating 
characteristics  of  the  systems  of  plans  involving  switching  between  normal  and 
tightened  inspection.  A  general  approach  to  determining  the  operating 
characteristic  of  sampling  schemes  using  Markov  chains  was  published  by  K.  S. 
Stephens  and  K.  E.  Larson  in  1967  [21].  They  applied  this  procedure  to 
investigate  the  behavior  of  the  105D  system  including  tightened,  normal  and 
reduced  inspection.  Insights  Into  the  behavior  of  the  sampling  scheme  as  a 
whole  were  provided.  The  article  Illustrated  the  benefits  of  incorporating 
the  system  in  the  following  ways:  (1)  affording  more  protection  for  the 
producer  as  well  as  the  consumer,  and  (2)  more  economies  of  sampling  effort. 
The  same  year  T.  L.  Burnett  [1]  also  developed  a  computer  program  for  the  IBM 
IBSYS  operating  system  to  compute  the  probability  of  acceptance  for  several 
sampling  schemes;  he  also  Imbeded  the  system  in  a  Markov  chain.  The 
calculations  are  then  reduced  to  matrix  manipulations  for  which  routines  were 
available;  his  research  Indicated  that  the  OC  curves  In  MIL-STD-105D  agree 
very  well  at  the  upper  end  of  the  curve  (region  of  smaller  values  of  fraction 
defective),  but  are  rather  pessimistic  on  the  lower  end  (regions  of  high 
fraction  defective  values)  when  compared  to  system  curves. 

A  computer  program  [20],  encompassing  the  modification  and  extensions  of 
the  previous  work  by  Stephens  and  Larson  [21]  and  Burnett  [1]  was  developed  in 
1977  by  J.  H.  Sheesley  at  the  General  Electric  Company,  and  was  run  on  a 
Honeywell  66/40  system;  the  program  evaluates  the  operating  characteristics 
for  the  Standard's  sampling  plans.  These  operating  characteristics  are 
obtained  through  the  use  of  Markov  chains  with  transition  probabilities  which 


describe  each  state  of  sampling  In  terms  of  the  state  which  was  occupied  at 
the  previous  sample. 


Curves  for  the  probability  of  acceptance  which  are  published  in  the 
standard  refer,  as  mentioned  earlier,  to  normal  Inspection  only;  the  curves 
produced  from  the  values  made  available  by  Sheesley's  program  justify  once 
again  the  necessity  and  importance  of  utilizing  the  switching  rules,  in  as 
much  as  the  use  of  only  the  normal  plans  sacrifices  protection  to  both  the 
producer  and  the  consumer  and,  furthermore,  also  requires  somewhat  more 
inspection  effort. 

Sheesley's  program  was  adopted  as  a  model  in  developing  the  program 
reported  herein,  while  the  approach  taken  and  procedures  for  scheme  evaluation 
were  patterned  after  the  work  of  Stephen  and  Larson  [21]. 

Some  of  the  noteworthy  publications  related  to  this  study  include  the 
works  of  H.  F.  Dodge  [4],  G.  J.  Keefe  [15],  W.  R.  Pabst,  Jr.  [19],  0.  A.  Cocca 
[2],  I.  D.  Hill  [14],  A.  J.  Duncan  [5],  B.  L.  Hansen  [13],  E.  L.  Grant  and  R. 
S.  Leavenworth  [7],  along  with  articles  by  A.  Hald  [9,  10],  R.  S.  Leavenworth 
and  R.  L.  Scheaffer  [17]  and  H.  C.  Hamaker  [11,  12]. 

Scope 

In  Chapter  II  the  methodology  and  the  approach  taken  in  developing  the 
computer  program  is  discussed.  Chapter  III  outlines  and  describes  the  program 
operation  and  capabilities  and  details  and  interprets  the  results  of  a  sample 


CHAPTER  2 
METHODOLOGY 


In  this  chapter,  the  approach  and  underlying  methodology  regarding  the 
design  and  modelling  of  the  techniques  for  evaluating  the  Standard  are 
presented.  Dr.  W.  R.  Pabst,  Jr.  had  this  to  say  while  considering  the  study 
of  the  operating  characteristics: 


.  .  .  they  follow  directly  from  the  statistical  parameters 
selected  and  are  uniquely  and  definitely  determined.  A 
sampling  scheme  is  the  whole  set  of  sampling  plans  and 
operations,  included  in  the  standard,  the  overall  strategy 
specifying  the  way  in  which  sampling  plans  are  to  be  used. 
[19:4] 


The  subject  of  calculating  and  evaluating  the  Standard  can  best  be 
described  by  the  following  comments  by  I.  D.  Hill: 


.  .  .  it  is  a  matter  of  art,  opinion,  aesthetic  sense  and 
compromise  as  well  as  of  science  and  mathematics  and  it  is 
this,  that  gives  the  subject  much  of  its  fascination.  .  .  . 
[14:31] 


At  this  stage,  proper  definition  of  the  term  "scheme"  in  conjunction  with  MIL- 
STD-105D  as  made  by  Stephens  and  Larson  [21]  is  used,  the  essence  of  which 
is  given  by  the  following  statement  attributed  to  I.  D.  Hill:  ".  .  .  an 
overall  strategy  specifying  the  way  In  which  sampling  plans  are  to  be  used" 
[14:31].  The  term  "system"  is  restricted  to  a  specific  set  of  normal, 
tightened,  reduced  inspection  plans  along  with  the  rules  for  switching  from 
one  to  another. 


Sampling  Scheme  Design 

The  generalized  form  of  the  switching  rule  criteria  are  shown  in  the  flow 
chart  in  Figure  2-1;  here  we  see  how  the  Standard  ties  together  sets  of  three 
sampling  plans,  each  at  a  different  level  of  severity,  into  a  unified 
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Figure  2-1.  FLOW  CHART  DEPICTING  THE  SWITCHING  RULE  CRITERIA 
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system.  Thus,  inspection  of  a  succession  of  lots  is  intended  to  move  among  a 
set  of  specified  tightened,  normal  and  reduced  sampling  plans  as  the  quality 
level  improves  or  degenerates.  The  Standard,  as  such,  allows  for  the 
application  of  the  normal  and  tightened  inspection  only,  while  reduced 
inspection  is  at  the  consumer  agency's  discretion.  Quality  levels  for  the 
producer  are  specified  in  terms  of  Acceptable  Quality  Level  (AQL),  while 


consumer  protection  is  afforded  by  the  switching  rules.  For  details  regarding 
the  specific  operation  of  the  switching  rules,  the  reader  is  referred  to  an 


article  by  A.  R.  Hahn  and  E.  G.  Schilling  [81  and  Paragraph  8.8  of  the 


Standard. 


Modelling  the  Operational  Behavior 


The  actual  behavior  of  the  process,  under  the  influence  of  the  sampling 


procedure,  can  be  very  dynamic  in  nature  due  to  shifts  in  the  severity  of  the 
inspection  levels.  To  evaluate  performance  of  the  process  under  such  condi¬ 
tions,  an  approach  similar  to  Stephens  and  Larson  [21]  is  adopted.  The 


probability  of  being  in  a  particular  state  depends  only  on  the  state  occupied 


at  the  previous  trial  or  inspection  and,  furthermore,  as  at  any  particular 


stage  of  a  sampling  plan,  the  state  of  the  scheme  is  completely  defined  by  one 


and  only  one  state. 


Thus,  we  see  that  the  sampling  plan  can  be  defined  in  terms  of  a  Markov 


chain.  We  can  then  obtain  the  steady  state  probability  of  being  in  each  state 


of  this  Markov  chain  in  terms  of  the  fraction  defective  ”p"  of  the  items  being 


inspected.  Hence,  for  each  stage  we  can  compute  such  quantities  as 


probability  of  accepting  the  lot,  number  inspected,  fraction  inspected  or 


other  characteristics  of  interest.  Once  these  state  probabilities  are 


obtained,  we  can  multiply  them  by,  say,  the  probability  of  acceptance  at  each 
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of  the  states,  yielding  the  overall  Probability  of  Acceptance  (PA);  similarly, 
the  ASN,  AOQ  and  the  AFI  also  may  be  obtained. 

Modelling  the  Evaluation  Procedure 
For  the  sake  of  completeness,  the  salient  features  regarding  the 
switching  rules  as  given  in  the  standard  [18,  5]  are  repeated  below: 

Normal  to  tightened:  When  normal  inspection  is  in  effect,  tightened 

inspection  shall  be  instituted  when  2  out  of  5 
consecutive  lots  or  batches  have  been  rejected 
on  original  inspection. 

Tightened  to  normal:  When  tightened  inspection  is  in  effect,  normal 

inspection  shall  be  instituted  when  5  consecu¬ 
tive  lots  or  batches  have  been  considered 
acceptable  on  original  Inspection. 


Normal  to  reduced:  When  normal  inspection  is  in  effect,  reduced 

inspection  may  be  instituted  when  on  normal 
inspection: 

(a)  The  preceding  10  lots  or  batches  (or  more 
as  indicated  by  the  note  to  Table  VIII) 
have  been  on  normal  inspection  and  none 
has  been  rejected  on  original  inspection; 
and 

(b)  The  total  number  of  defectives  (or 


defects)  in  the  samples  from  the 
preceding  10  lots  or  batches  (or  such 
numbers  as  was  used  in  condition  a  above) 


is  equal  to  or  less  than  the  applicable 
[limit!  number,  "M",  given  in  Table  VIII 
[of  the  Standard!.  If  double  or  multiple 
se  pling  is  in  use,  all  samples  inspected 
should  be  Included,  not  "first"  samples 
only;  and 

(c)  Production  is  at  a  steady  rate;  and 

(d)  Reduced  inspection  is  considered 
desirable  by  the  responsible  authority. 

Reduced  to  normal:  Inspection  shall  revert  back  to  normal 

inspection  if  any  of  the  following  conditions 

occur: 

(a)  A  lot  or  batch  us  rejected; 
or 

(b)  A  lot  is  marginal,  i.e.,  the  number  of 
defectives  (or  defects)  falls  between  the 
applicable  acceptance  and  the  rejection 
numbers;  o£ 

(c)  Production  becomes  irregular  or  delayed; 
or 

(d)  Other  conditions  warrant  that  normal 
inspection  be  reinstated. 

The  reader  is  referred  to  Figure  2-1  for  an  overview  of  the  application 


of  the  switching  rules. 


Development  of  a  solution  method  requires  the  following  steps: 

1.  Drawing  of  a  state  transition  diagram,  shown  in  Figure  2-2,  con¬ 
taining  a  state  for  each  stage  required  for  the  operation  of  the  sampling 
scheme  system.  The  interconnecting  arrows  represent  all  permitted  transitions 
from  one  state  to  another  for  all  stages.  A  state  is  defined  by  ,  i  * 

1,...,  20. 

2.  Assigning  to  each  transition  an  appropriate  probability.  P^j  is  the 
one  step  transition  probability  of  going  from  state  to  Sj;  for  any  two 
states  not  connected  by  an  arrow  P^-  =  0. 

3.  Obtaining  the  Probability  Transition  Matrix,  Figure  2-3,  from  the 
state  transition  diagram.  Figure  2-2. 

4.  Solving  the  above  matrix  so  that  the  steady  state  probabilities  of 
being  on  tightened,  normal,  and  reduced  inspection  can  be  found. 

5.  Finally,  multiplying  appropriate  combinations  of  the  probabilities 
obtained  in  Step  4  by  the  characteristics  of  each  state  in  order  to  obtain 
such  characteristics  as  ASN,  AFI,  and  AOQ. 

In  Step  1,  the  state  transition  diagram  presents  a  flow  diagram  of  the 
operation  a  system.  Five  states,  SI,  S2,  S3,  S4,  S5,  are  required  for 
tightened  inspection  (if  five  consecutive  lots  are  accepted,  then  normal 
inspection  is  instituted).  Next,  14  states,  S6,  S7,...,  S19,  handle  the 
various  cases  for  normal  inspection.  Ten  states  are  required  to  represent  10 
consecutive  lots  being  accepted  and  thus  heralding  reduced  inspection;  four 
states  handle  the  situation  where  two  out  of  five  consecutive  lots  have  been 
rejected  and  inspection  is  switched  to  the  tightened  level;  finally,  one 
state,  S20,  represents  the  reduced  inspection  level,  where  a  single  rejected 
lot  reverts  inspection  back  to  the  normal  level. 
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Figure  2-3.  TRANSITION  PROBABILITY  MATRIX  FOR  MIL-STD-105D 


In  Step  2,  appropriate  probabilities  are  assigned  to  the  various 
transitions.  All  the  mathematical  work  at  this  step  and  also  throughout  this 
study  is  based  on  the  binomial  probability  function.  The  probabilities  are  as 


I' 


■s 


i 


* 


j 


fc; 


S 


a 


follows: 


T  =  Probability  of  being  accepted  on  tightened  inspection. 

N  =  Probability  of  being  accepted  on  normal  inspection. 

R  =  Probability  of  being  accepted  on  reduced  inspection. 

Q  =  Probability  of  no  more  than  M  defectives  in  the  last  10  lots 

when  on  normal  inspection. 

n-i 


T  * 


r  J  „r 

E  (,•)  p  q 

i=0 


N  = 


E  (?)  P1  q 
i=0 


n-i 


R  = 


v  in\  o1 
E  ( i )  p  q 
i=0 


M 


where  in  general: 


1  (  «  )  P  q 
i=0  1 


c 

n 

M 


acceptance  number 
sample  size 

limit  number  for  reduced  inspection  ([lfil  Table 
VIII) 


Step  (b)  of  the  criteria  for  switching  from  reduced  to  normal  inspection 
requires  that  no  more  than  a  certain  number  of  defectives  or  (defects),  "on" 
be  found  in  10  consecutive  lots  inspected.  The  probability  of  passing  this 
limit  number  criterion  is  approximated  by  the  binomial  probability  of  finding 
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a  count  less  than  or  equal  to  the  prescribed  number  of  defectives  (M)  in  a 
sample  10  times  the  reduced  inspection  sample  size  (or  the  ASN,  in  the  case  of 
double  sampling  plans).  This  approximation  is  attributed  to  T.  L.  Burnett 
[1].  This  situation  is  dealt  with  in  more  detail  by  Stephen  and  Larson  [211, 
but  one  cannot  arrive  at  a  closed  form  analytical  expression  from  their  work 
thus  requiring  estimation  by  iteration.  J.  H.  Sheesley  confirms  the 
indeterminate  form  of  this  event  and  states:  ".  .  .  in  any  event  the  Burnett 
approximation  should  provide  reasonable  accuracy"  [20:ii1.  A  graphical 
comparison  of  OC  curves  for  the  sample  plan  for  code  letter  L  at  1%  AQL, 
obtained  by  both  T.  L.  Burnett  [1]  and  by  Stephens  and  Larson  [211,  is 
presented  in  Appendix  A;  readers  are  referred  to  the  appropriate  papers  for 
further  details. 

Step  3  describes  the  transition  probability  matrix  for  the  sampling 
system  under  investigation.  The  construction  is  such  that  the  first  five  rows 
indicate  the  event  of  the  process  being  on  tightened  inspection,  where  five 
accepted  lots  in  a  row  are  needed  to  switch  to  normal  inspection:  the  next 
four  rows  represent  states  in  which  normal  inspection  is  in  effect  and  a 
rejection  followed  by  up  to  three  acceptances  has  occurred:  when  the  rejection 
is  followed  by  four  acceptances  the  criterion  of  two  out  of  five  rejections 
(see  Figure  2-1)  causing  a  switch  to  tightened  inspection  no  longer  applies. 

Th  next  10  rows  represent  the  possibility  of  up  to  nine  consecutive  lots  being 
accepted  on  normal  inspection.  At  row  19  the  possibility  of  passing  the  limit 
number  criterion,  as  discussed  previously,  is  represented;  in  case  the  lot  is 
accepted  but  the  limit  number  requirement  is  not  passed,  then  the  procedure  is 
designed  such  that  the  criterion  of  being  on  normal  inspection  is  switched 
back  to  the  first  lot  in  the  sequence;  this  procedure  keeps  in  check  the 
possibility  of  having  a  slightly  liberal  estimate  for  the  probability  of  being 


on  reduced  inspection.  Row  20  shows  all  possible  occurrences  while  on  reduced 
inspection. 

Finally,  once  the  state  probabilities  are  obtained  then,  using  the  method 
detailed  above,  the  overall  system  steady  state  probabilities  can  be 
obtained.  The  next  chapter  describes  and  details  the  computer  program 
developed  to  evaluate  the  above  steps. 


rHAPTFR  X 

PROGRAM  DESCRIPTION 


The  program  Is  written  in  structured  FORTRAN  IV  and  was  coded  and  run  on 
a  VAX  11/750  machine,  working  under  the  VMS  operating  system  and  utilizing  an 
ADM-3A  Interactive  display  terminal.  The  program  is  completely  interactive 
and  may  be  described  as  allowing  users  to  input  information  in  the  same  order 
as  they  would  use  to  select  a  system  of  plans  from  the  Standard. 


Programming  Style 


In  developing  this  user-friendly  program,  all  efforts  were  made  to  adhere 
to  a  number  of  rules  and  conditions  applicable  to  writing  of  interactive 
routines  as  suggested  by  Gaines  and  Facey  [6].  Users  are  expected  to  follow 
the  program  model  step  by  step  and  answer  questions  to  make  choices  as  to 
their  information  needs,  parameter  usage,  etc.  Thus,  a  novice  operator  may 
follow  the  structure  of  the  program  and  have  no  trouble  in  gettin  through  the 
maze. 

The  program  is  called  MILSTD.  The  following  list  explains  the  function 
of  the  different  elements  or  subprograms  of  MILSTD: 


LIST  OF  PROGRAM  ELEMENTS  AND  THEIR  FUNCTION 


SUBPROGRAM 


FUNCTION 


Identifies  the  code  letter  associated 
with  any  lot  size  for  a  particular 
inspection  level  (Table  I  [18]). 


1 

k  k  f!  vv*.  t  r  ww; 

p 

SUBPROGRAM 

FUNCTION 

i  * 

SS 

Find  the  sample  size(s)  for  either  the 

h: 


l 


I 


INDEX 


VALUES 


OC 


PLOT  OC 


single  or  the  double  sampling  plan  for 
the  code  letter  obtained  from  CODE 
(Tables  II-A,  B,  C  and  III-A,  B,  C 
[18]). 

Locates  the  particular  cell  in  the 
tables  of  the  Standard  for  the 
designated  AQL  and  code  letter. 

Obtains  the  acceptance  and  rejection 
numbers  for  either  the  single  or  double 
sampling  plan  for  the  cell  located  by 
INDEX  (TABLES  II-A,  B,  C  AND  III-A,  B, 

C  [18]). 

Evaluates  the  sampling  plan  and  obtains 
the  Probability  of  Acceptance  (PA), 
average  sample  number  (ASN),  Average 
Outgoing  Quality  (AOQ)  and  the  Average 
Fraction  Inspection  (AFI)  for  each 
value  of  fraction  defective,  p, 
entered. 

Provides  plots  of  the  OC,  AOQ,  ASN.  and 
AFI  curves  for  the  system. 
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Program  Operation 

The  operation  starts  with  the  loading  of  MILSTD:  the  user  supplied 
requirement  (inputs)  are  as  follows: 

1.  Inspection  level:  The  desired  inspection  level,  "SI,"  "S2," 

"S3,"  "S4“  for  special  and  "1,"  "2"  or  "3"  for 
general  inspection. 


2.  Lot  size: 

I 


The  actual  number  of  items  in  the  lot 
being  sampled. 


3.  AQL:  The  designated  AQL  is  entered  in  percent 

and  must  be  one  of  the  prescribed  values  given 
in  the  Standard,  otherwise  strange  results  may 
be  obtained. 

4.  Sampling  plan:  The  user  is  allowed  the  option  to  select 

either  a  single  or  double  sampling  plan,  the 
choice  can  be  indicated  by  entering  "S"  or  "D" 
for  the  appropriate  plan. 

5.  Evaluation  request:  The  program  asks  the  user  if  evaluation  of 

the  plan  presented  so  far  is  required.  If  the 
decision  is  yes,  enter  1;  if  not,  enter  2. 

6.  Table  or  Graph  format:  The  program  requires  the  user  to  choose 

the  style  of  output. 
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The  following  two  inputs  are  required  only  if  the  selection  at  input  6 
was  for  a  tabulated  format: 

7.  Number  of  fraction  defective  values  for  which  the  evaluations  are  to 
be  made. 

8.  Values  of  the  fraction  defective. 

Once  the  user  has  provided  the  values  or  the  decision  criteria,  as  the 
case  may  be,  the  program  returns  a  completely  evaluated  MIL-STD-105D  sampling 
system  in  either  a  tabulated  or  graphical  form.  Data  for  inputs  1  and  4  must 
be  placed  in  apostrophy  marks.  Input  5  is  required  in  the  sense  that,  by 
then,  the  user  has  a  summary  of  the  sampling  plan  for  the  normal,  tightened 
and  reduced  inspection  levels,  which  consists  of  the  sample  size,  acceptance 
number,  and  rejection  number  for  each.  Alternatively,  the  program  may  display 
the  message:  "USE  100  PERCENT  INSPECTION  AS  SAMPLE  SIZE  EQUALS  OR  EXCEEDS  LOT 
SIZE."  In  the  latter  case,  scheme  evaluation  does  not  occur  and  the  program 
terminates.  Another  situation,  which  may  occur  when  double  sampling  plans  are 
desired,  is  summarized  in  the  following  self-explanatory  message:  "SAMPLE 
SIZE  NOT  AVAILABLE  AT  ANY  INSPECTION  LEVEL,  CORRESPONDING  SINGLE  SAMPLING  PLAN 
SHOULO  BE  USED."  The  program  then  displays  the  appropriate  single  sampling 
plan.  On  receiving  this,  the  user  can  exercise  the  option  of  either 
evaluating  the  single  sampling  system  under  investigation  or  stopping  the 
run.  This  technique  of  offering  a  number  of  simple  choices  and  explicit 
instructions  for  carrying  out  commands  is  known  as  a  "structured  walk¬ 
through."  Figure  3-1  illustrates  the  program  structure. 

The  program's  computational  behavior  is  illustrated  by  flow  charts  in 
Figure  3-2.  The  sequential  operation  of  the  different  processes  becomes  clear 


Figure  3-1.  PROGRAM  ORGANIZATION  CHART 
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START 


Figure  3-2.  FLOW  CHART 


as  one  relates  this  flow  chart  to  a  sample  run  for  a  typical  sampling  plan. 

The  sample  run  is  discussed  in  the  following  section. 

Program  Capability  (or  the  Sample  Run) 

The  sample  run  looks  at  a  double  sampling  system  where  the  user  is 
working  with  a  lot  size  of  1000,  Inspection  Level  "II,"  and  a  designated  AOL 
of  2.5%.  The  run  begins  with  the  user  entering  the  appropriate  selections. 

The  plan  is  summarized  for  acceptance  (AC)  and  rejection  (RE)  numbers,  as  well 
as  sample  size  for  both  the  stages  and  all  three  sampling  plans  of  the  system 
(see  Appendix  C  for  an  illustration). 

The  user  is  then  asked  to  either  stop  the  run  or  continue  with  the 
evaluation  as  the  particular  situation  may  warrant.  Once  a  choice  is  made  to 
progress  further  with  the  run,  the  user  is  at  liberty  to  select  either  a 
tabulated  format  for  one  or  more  specific  values  of  fraction  defective  "p,"  or 
a  graph  format  whereby  curves  are  plotted  for  the  0C,  ASN,  AOQ  and  AFI  for  the 
full  range  of  "p".  The  output  of  this  sample  run,  for  both  formats,  are 
presented  in  Appendix  C.  The  computer  printer  must  be  capable  of  being  set  to 
compressed  printing  for  a  clear  output  of  the  graphs. 

Analysis  and  Interpretations 

The  study  of  the  tabulated  and  plotted  values  for  the  operating 
characteristics  of  the  Standard  based  on  full  utilization  of  the  switching 
rules  leads  to  some  interesting  observations,  which  are  detailed  as  follows: 

1.  There  is  a  marked  improvement  in  the  0C  curve  in  comparison  to  that 
given  for  the  normal  inspection  only  brought  about  by  using  the  switching 
rules.  From  Figure  3-3  it  is  seen  that  the  risks  of  accepting  product  of 
quality  worse  than  the  desired  AQL  are  drastically  reduced. 
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Figure  3-3.  OPERATING  CHARACTERISTIC  CURVES,  CODE  J,  2.5%  AQL 
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2.  When  rectifying  inspection  is  employed,  the  Average  Outgoing  Quality 
Limit  for  the  scheme  is  much  improved  over  the  levels  reported  for  the  normal 
inspection  plan  only  as  shown  in  Figure  3-4.  The  A00  of  the  product  is 
important  to  the  consumer  because  it  represents  the  average  quality  ultimately 
received.  The  rectification  features  of  the  scheme  ensure  an  AOQ  as  good  as, 
and  more  often  better  than,  the  incoming  quality. 

3.  In  the  case  of  the  double  sampling  system,  different  sample  sizes  are 
involved  for  code  letter  J  at  2.5%  AQL  (or  for  that  matter,  any  other 
combination),  the  sample  size  for  the  scheme  can  then  only  be  represented  as 
an  expected  value.  This  is  the  ASN  for  the  scheme.  Figure  3-5  shows  the 
reduction  in  the  ASN  in  the  region  where  quality  is  good,  since  there  is  a 
possibility  of  going  to  reduce  inspection.  At  regions  where  the  incoming 
fraction  defective  is  somewhat  higher,  the  second  sample  is  drawn  more  often 
and  hence  the  ASN  increases,  indicating  the  institution  of  tightened 
inspection. 

4.  The  AFI  curves  are  given  as  an  indication  and  guide  to  determine  the 
inspection  effort  involved  when  rectifying  inspection  is  used.  As  can  readily 
be  seen  in  Figure  3-6,  there  is  a  reduction  in  the  AFI  for  good  quality- 
regions  and  a  better  check  is  imposed  than  the  normal  inspection  plan  for 
higher  values  of  fraction  defective. 

Thus,  it  can  be  seen  that  the  producer  of  good  quality  products  can  seek 
the  reward  of  switching  to,  and  remaining  on,  reduced  inspection,  in  terms  of 
lower  sample  sizes  and  less  inspection  costs,  yet  the  protection  afforded  to 
the  consumer  remains  the  same,  in  as  much  as  the  reduced  inspection  is 
instituted  and  significant  at  only  the  lowest  level  of  percentage  defective. 


CHAPTER  4 
CONCLUSION 


This  project  had  two  purposes.  One  was  to  introduce  the  MIL-STD-105D 
into  the  computer  as  a  user  friendly,  interactive  package,  considering  single 
and  double  sampling  plans  only.  The  second  was  to  provide  analytical  output 
of  characteristics  and  effectiveness  measures  of  the  systems  of  normal, 
tightened  and  reduced  inspection  plans  in  either  tabular  or  graphical  form. 


Summary 


In  Chapter  2  it  was  shown  that  the  different  sampling  plans  of  the 
Standard  formed  a  sampling  system,  and  how  this  overall  system,  including  the 
application  of  the  switching  rules,  could  be  described  by  a  Markov  chain. 
Utilization  of  the  Markovian  property  was  made  in  the  sense  that  the  past 
history  was  completely  summarized  in  the  specification  of  the  current  state, 
thereby  a  matrix  of  transition  probabilities  could  be  derived.  Thus,  an 
approach  was  taken  to  solve  this  Markov  problem  which  gave  the  steady  state 
probabilities  of  being  on  normal,  tightened  or  reduced  inspection.  These 
steady  state  probabilities  allowed  the  formulation  of  the  system 
characteristics,  hence  a  methodology  to  evaluate  the  system  was  developed. 

Finally,  the  software  was  described  in  Chapter  3.  The  user  should  be 
able  to  understand  the  program  features  and  capabilities  as  well  as  how  the 
various  routines  of  the  program  integrate. 
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TABLE  A-l. 

CODE  LETTER:  L  AQL:  1  % 

P 

Stephen  and 

Larson's  Method 

Burnett's 

approximation 

0.002 

1.00 

1.00 

0.004 

1.00 

0.99 

0.006 

1.00 

0.99 

0.008 

0.99 

0.99 

0.98 

0.98 

0.012 

0.95 

0.95 

0.014 

0.88 

0.89 

0.016 

0.75 

0.76 

0.018 

0.55 

0.60 

0.020 

0.45 

0.46 

0.030 

0.14 

0.14 

0.040 

0.04 

0.04 

0.050 

0.01 

0.01 
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Figure  A-l.  OC  CURVE  FOR  SAMPLE  PLAN  "L"  AT  AQL  =1.0% 
STEPHENS  AND  LARSON  METHOD 
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Figure  A-2.  OC  CURVE  FOR  SAMPLE  PLAN  "L"  AT  AQL  =  1.0% 
BURNETT'S  APPROACH 


APPENDIX  B 

LIST  OF  KEY  VARIABLES  USED  IN 
THE  PROGRAM  AND  THEIR  FUNCTION 
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VARIABLE 


DEFINITION 


LOCATION 


AOL  Acceptance  Quality  Level  Main  program. 

Subroutine  LNM 


COL 

FNR 

FRR 

FS 

FTS 

I 

J 

L 

LNC 

NFA 

NNNN 

NQ 

NTN 

NZ 

P 

PA 


Code  letter 

Rejection  number  for  normal 
inspection  at  the  first  sampling 
stage 

Rejection  number  for  reduced 
inspection  at  first  sampling 
stage 

First  stage  sample  size  for 
normal  or  tightened  inspection 

Rejection  number  for  tightened 
inspection  at  the  fist  sampling 
stage 

Inspection  level 
Number  of  stages 
Lot  size 

Limit  number  for  reduced 
inspection  "M" 

Acceptance  number  for  normal 
inspection  at  the  first  sampling 
stage 

Lot  size 

Sample  size  at  normal  or 
tightened  inspection  for  single 
sampling  plan 

Acceptance  number  for  normal 
inspection  at  the  second 
sampling  stage 

Ten  times  the  average  sample 
size  at  normal  inspection,  used 
to  calculate  the  limit  number  "M" 

Value  of  fraction  defective 


Subroutine  code 
Subroutine  values 

Subroutine  values 

Subroutine  SS 
Subroutine  values 

Main  program 

Subroutine  OC 

Main  program 

Subroutine  OC, 
Subroutine  LNM 

Subroutine  values 

Subroutine  OC 
Subroutine  SS 

Subroutine  values 

Subroutine  OC, 
Subroutine  LNM 

Subroutine  OC 


Probability  of  Acceptance 


Subroutine  OC 


VARIABLE 


DEFINITION 


LOCATION 


0 

Probability  of  passing  the  limit 
number  criteria 

Subroutine  OC 

RFA 

Acceptance  number  for  reduced 
inspection  at  the  second 
sampling  stage 

Subroutine  values 

RFS 

Sample  size  for  reduced  inspec¬ 
tion  at  second  sampling  stage 

Subroutine  SS 

RN 

Sample  size  at  reduced  inspec¬ 
tion  for  single  sampling  plan 

Subroutine  SS 

RRI 

Probability  of  being  accepted  on 
reduced  inspection 

Subroutine  OC 

APPENDIX  C 

OUTPUT  OF  BOTH  TABLE  AND  GRAPH 
FORMATS  FOR  THE  SAMPLE  RUN 


ENTER  INSPECTION  LEVEL  IN  QUOTATION  MARKS 
E.  G.  »  SPECIAL  :  'SI'.  'S2 ' .  'S3'.  'S4' 

GENERAL  :  '1  '.  '2'.  '3' 


ENTER  LOT  SIZE: 

1000 

ENTER  AQL  IN  PERCENT 
2.  5 

DO  YOU  WANT  SINGLE ('S')  OR  DOUBLE ('D') 
SAMPLING  PLANS  I (NOTE: ENTER  S  OR  D  IN  QUOTES). 
'D' 


THESE  PLANS  ARE  : - 
*:  NORMAL: :::::::::: 


TIGHTENED:  :::::::::::: REDUCED: : * 


AC  1  = 
RE  1  = 
Nl  = 

2 

5 

50 

AC  1  = 

RE  1  = 

Nl  = 

1 

4 

50 

AC  1= 

RE  1  = 

Nl  = 

0 

4 

20 

AC  2= 

6 

AC  2— 

4 

AC  2= 

3 

RE  2= 

7 

RE  2- 

5 

RE  2- 

6 

N2- 

50 

N2= 

50 

N2= 

20 

DO  YOU  WANT  SCHEME  EVALUATION. . ? 

IF  YES  ENTER . 1 

IF  NO  ENTER . 2 

1 

DO  YOU  WANT  A  TABLE  OR  A  GRAPH  FORMAT  7 
FOR  GRAPH. ...  ENTER:  1 
FOR  TABLE _ ENTER:  2 

2 

SPECIFY  THE  NUMBER  OF  FRACTION  DEFECTIVE  VALUES. 
4 

ENTER  THE  FRACTION  DEFECTIVE  VALUES (S). 

(PUT  A  COMMA  BETWEEN  VALUES.  > . 

•  Ol#  •  OS#  •  OS#  •  10 


SCHEME  OPERATING  CHARACTERISTICS 


P:  :  ::::: 

:  :  :  P  ( A ) : :  :  : : 

:  :  ASN:  :  :  : 

:  :  AOQ: :::::: 

:  :  : AFI 

0.  OlO 

1.0000 

23.  82 

0.  0098 

0.02 

0.  050 

0. 5242 

73.  17 

0.  0244 

0.  51 

0.  080 

0. 1316 

67.  13 

0.  0098 

0.  88 

0.  100 
FORTRAN 
* 

0. 0472 
STOP 

60.  83 

0.  0044 

0.  96 

O’v'v  'j 


V'l’.  *.«i* 


CV.CU.NiC>.  (.(.I. 


^.jV^VWV*  *  A  AAA1  -.*  •  *  o  ■/ *l<  •  n 


'*  W ■-’  Oa*"  *.  '* 


MMM 
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MILSTD. FOR 

THIS  AN  INTERACTIVE  PROGRAM  TO  EVALUATE 
THE  MIL-STD— 10SD  SAMPLING  SCHEME. 

PROGRAMMED  BY:  AZMAT  H. SIDDIQI 
INDUSTRIAL  AND  SYSTEMS  ENGINEERING  DEPT. 
UNIVERSITY  OF  FLORIDA. GAINESVILLE. 
FLORIDA  32611 


PROGRAM  MAIN 
CHARACTER  I *2. CDL.  SP 

INTEGER  FS.  FA.  SA.  AR,  SR.  AC.  AT.  AE.  SC.  ST.  Z.  A.  B 
INTEGER  NN,  NFS.  NFA.  NAE.  NSA.  NAR.  NAC.  NAT.  TN. TFS. TFA.  TAE.  TSA 
INTEGER  TAR.  TAC.  TAT.  RN.  RFA. RAE.  RAR.  RAC.  RAT. FRR.  FTR.  FNR.  SNR 
«  .  STR.  SRR.  RFS.  SI.  S2,  S3.  S4.  S 

COMMON/SCL/L.  SP.  S.  I.  NO.  T.  R.  CDL.  FS.  D.  J.  K.  AQL.  M.  FA.  SA. 

•  AR. SR. AC. AT. AE. SC. ST.  Z.  NTN.  NFS.  NFA.  NAE. NSA. NAR.  NAC.  NAT. 

*  TN.  TFS. TFA.  TAE.  TSA.  TAR.  TAC.  TAT.  RN.  RFS. RFA. RAE.  RAR.  RAC.  RAT 
»  . FRR. FTR.  FNR.  SNR.  STR.  SRR.  NNNN 


THIS  SECTION  ASKS  THE  USER  TO  ENTER  VALUES  FOR 
THE  NECESSARY  PARAMETERS. 


PRINT  *.  'ENTER  INSPECTION  LEVEL  IN  QUOTATION  MARKS  ' 
PRINT  *.  'E.  C.  I  SPECIAL  :  '  'SI  *  *.  "S2  ".  "S3".  "S4  "  ' 
PRINT*.'  GENERAL  :"1".  "2".  "3"  ' 

READ  *.  I 

PRINT  *. 'ENTER  LOT  SIZE: ' 

READ  *.  L 

PRINT  *. 'ENTER  AQL  IN  PERCENT' 

READ  *. AQL 

PRINT  *.  'DO  YOU  WANT  SINGLE!  "S"  >  OR  DOUBLE*  "D">  * 
PRINT  *.  'SAMPLING  PLANS  >  (NOTE: ENTER  S  OR  D  IN  QUOTES).  ' 
READ  *. SP 
PRINT*.  X  ' 

PRINT  «.  'THESE  PLANS  ARE  ' 

CALL  CODE 
CALL  SS*» 

CALL  INDEX 
CALL  VALUES 
S-O 

IF(RFA.LT. 100) THEN 
FRR -RFA* 1 
END  IF' 

IFtRACILT. 100) THEN 
RFA— RAC 
FRR-RFA+3 
ENDIF 

IF (RAT. LT. 100)THEN 
RFA— RAT 

FRR-RFA+2 
ENDIF  * 

IF(RAE.LT. 100) THEN 
RFA— RAE  . 

FRR-RFA+3 

ENDIF 

IF(RAR. LT. 100>THEN 

RFA-RAR 

FRR-RFA+4 

ENDIF 

IF(TFA.  LT.  100) THEN 

FTR-TFA+1 

ENDIF 

IF* TAT. LT. 1 OO ) THEN 

TFA-TAT 

FTR-TFA+2 

ENDIF 

IF(TAC.LT. 100) THEN 

TFA-TAC 

FTR-TFA+3 

ENDIF 

IF ( TAR. LT. 100) THEN 
TFA- TAR 
FTR-TFA+4 
ENDIF 

IF (TAE. LT. 100)THEN 


oonnonnon 


TFA-TAE 
FTR-TFA-O 
END  IF 

IF1NFA.LT. 100 ) THEN 

FNR-NFA«-1 

END  IF 

IF1NAT. LT. 100) THEN 
NFA-NAT 
FNR-NFA+2 
END  IF 

IF1NAC.LT. 100) THEN 
NFA-NAC 
FNR-NFA+3 
END  IF 

IF1NAR.LT.  1 00  >  THEN 
NFA-NAR 
FNR-NAR+4 
END  IF 

IFINAE.LT.  100 ) THEN 
NFA-NAE 
FNR-NFA+3 
END  IF 

IF1NSA.LT. 100) THEN 
NTN-NSA 
SNR-NTN+1 
END  IF 

IF<  TSA. LT. 100) THEN 

TN-TSA 

STR-TN-M 

END  IF 

IF(SC.LT. 100) THEN 

Z-SC 

SRR-Z+3 

END  IF 

IF1ST.LT.  100) THEN 

Z-ST 

SRR-Z+2 

END  IF 

IF1SR.LT. 100) THEN 
Z-SR 
SRR-Z+4 
ENDIF 

IF (RFS.  CE.  L.  OR.  FS.  GE.  L.  OR. NO. CE.  L.  OR.  RN.  CE.  L ) THEN 

this' section'  contains'  the' format' statements . :  : 

NEEDED  TO  SET  UP  TABLES  FOR  THE  PLAN  STATISTICS  : 

AT  THE  DIFFRENT  INSPECTION  LEVELS. 

THE  USER- IS  ASKED  IF  SCHEME  EVALUATION  IS 
DESIRED.  JF  SO  THEN  SUBROUTINE  OC  IS  INVOKED.  :: 

PRINT  *.  'USE  lOO  PERCENT  INSPECTION.  ' 

PRINT  *.  'AS  SAMPLE  SIZE  EQUALS  OR  EXCEEDS  LOT  SIZE. ' 

S-l 

ENDIF 

IF1S.EQ.  l)CO  TO  60 

IFUSP.EQ.  'S').  OR.  (SP.EQ.  '0'.  AND.  CDL.NE.  'A'MTHEN 

PRINT  *.  NORMAL: :::::::::::: TIGHTENED: :::::::::::: REDUCED: : * 

print  *,  ' - - - 

WRITE (3.  10)NFA. TFA.  RFA 

10  FORMATOX.  'AC  1»'»  13.  10X.  'AC  1-M7.  I0X.  'AC  1-'.  17) 

WRITE13.  201FNR.  FTR.  FRR 

20  FORMATOX.  'RE  1»'.  13.  10X,  'RE  1-'.  17.  10X.  'RE  1-'.  17) 

ENDIF 

IF1SP.EQ. 'D'. AND. CDL.NE. 'A'lTHEN 
WRITE (3. 23)FS.  FS.  RFS 

23  FORMATOX.  'Nl»'»  13.  10X.  'Nl«'.  19,  10X,  'Nl-'.  19) 

PRINT  «.  ' - 

ENDIF 

IF1SP.E0. 'D'. AND. CDL.NE. 'A'lTHEN 
WR I TE 1  3.  30 )  NTN.  TN.  Z 

30  FORMATOX.  'AC  2-'.  13.  10X.  'AC  2-'.  17.  IOX.  'AC  2-'.  17) 

WRITEO,  40) SNR.  STR.  SRR 

40  FORMATOX.  'RE  2-'.  13.  IOX.  'RE  2-'.  17.  IOX.  'RE  2-'.  17) 

WRITEO.  43)FS.  FS.  RFS 

43  FORMATOX.  'N2-'.  IS.  IOX.  'N2-'.  19.  IOX.  'N2-'.  19) 

ENDIF 

PRINT  *.  ' - 

IF1SP.E0.  'S') THEN 
WRITEO.  30>N0.  NO.  RN 

30  FORMATOX.  'N-'.  16.  IOX.  'N-'.  110.  IOX.  'N«'.  110) 

PRINT  *. ' - 


*■*  »,r*T)T» 
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IfCCSP.EQ.  'D'.AND.CDL.EQ.  'A'.AND.  J.CE.  19>.0R.  CSP.EO.  'D'.AND.CDL. 

*  PRINT  if^SAMPLl‘siZETNOT  AVAILABLE  AT  ANY  INSPECTION  LEVEL*  * 
PRINT  *.  'CORRESPONDING  SINGLE  SAMPLING  PLAN  WILL  BE  USED: 

PRINT  *.  '  ' 

SP«'S' 

GO  TO  1 

PRINT*. 'DO  YOU  WANT  SCHEME  EVALUATION..?' 

PRINT  #.  'IF  YES  ENTER . 1' 

PRINT  *.  'IF  NO  ENTER . 2' 

READ  *.  AZ 

IF<AZ. EQ. 2)  CO  TO  60 
IF (AZ. EO. 1>  GO  TO  38 
38  NNNN-L 

IF  C SP. EQ. 'S')  JT-1 
IF ( SP  EQ  ji‘*2 

CALL  OCCNNNN. FS.  RFS.  RFA.  SRR.  FRR.  NO.  RN.  TFA.  TN.  FTR.  STR.  NFA.  NTN 

*  . FNR. SNR. Z. JT. AOL> 

60  STOP 

END 

SUBROUTINE  SS 
C 

C***************************************************** 

C**  THIS  SUBROUT I NE  IF I NDS  THE  SAMPLE  SIZES  FOR  ** 

C**  BOTH  THE  DOUBLE  AND  THE  SINGLE  SAMPLING  PLANS  ** 

C**  AT  NORMAL. TIGHTENED  AND  REDUCED  INSPECTION  ** 

C**  LEVELS.  ** 

C***************************************************** 

c 

INTEGER'S.1  fa!  I/l'aR.  SR.  AC.  AT.  AE.  SC.  ST.  Z.  NN.  NFS.  NFA.  NAE.  NSA 
INTEGER  NAC. NAT.  TN.  TFS.  TFA.  TAE.  TSA.  TAR.  TAC. TAT.  RN. RFA. RAE 
INTEGER  RAR.RAC.  RAT.  NAR.FRR.RFS.  FTR.  FNR.SNR. STR.  SRR 
COMMON/SCL/L.  SP.  S.  I.  NO.  T.  R.  CDL.  FS.  D.  J.  K. AQL. M. FA. SA. AR 
«  .SR.AC.AT.AE.SC.  ST.  Z.  NTN.  NFS.  NFA.  NAE.  NSA.  NAR. 

*  NAC.  NAT.  TN.  TFS.  TFA.  TAE.  TSA.  TAR,  TAC.  TAT. RN. RFS. RFA 
C  . RAE.  RAR. RAC.  RAT,  FRR.  FTR,  FNR.  SNR,  STR.  SRR 

IF<SP.EG. 'S'LTHEN 

IFCCDL. EQ.  ‘ A ' >  THEN 

NQ?2 

EN&IF 

IF < CDL. EQ. 'B'lTHEN 

NQ-J 


END  IF 

IFCCDL. EQ. 'C ' )  THEN 

NO-5 

END  IF 

IFCCDL.EQ.  'D')  THEN 

NQ-B 

END  IF 

IFCCDL.EQ.  'E ' ) THEN 
NO- 13 
END  IF 

IFCCDL.EQ. 'F'JTHEN 
NO— 20 
END  IF 

IFtCDL. EQ. 'G ' ) THEN 

NQ-32 

END  IF 

IFCCDL.EQ. 'H'»  THEN 
NO- 50 
END  IF 

IFCCDL.EQ.  ' J ' > THEN 
NO— 80 
END  IF 

IFCCDL.EQ. 'K ' ) THEN 

NQ-125 

END  IF 

IFCCDL.EQ. 'L ' ) THEN 
NO— 200 
END  IF 

IFCCDL.EQ.  'M ' 1  THEN 
NQ-3 13 
END  IF 

IFCCDL.EQ.  'N ' ) THEN 
NO* 500 
END  IF 

IFCCDL.EQ.  'P ' 1  THEN 

NQ-800 

END  IF 

IFCCDL.EQ.  '0  ) THEN 


ENDIF 


Nu* 1230 
ENDIF 

IF(CDL.EO. 'R  *  > THEN 

NQ-2000 

ENDIF 

IF(CDL.EQ.  'S' ) THEN 

NQ-3130 

ENDIF 

IF(CDL.EQ.  'A'.OR.CDL. EQ.  'B  '. OR. CDL. EQ.  'C'JTHEN 

RN-2 

ENDIF 

IF(CDL. EQ.  'D'l  THEN 

RN-3 

ENDIF 

IF (CDL. EQ.  'E'»  THEN 

RN-5 

ENDIF 

IF (CDL. EQ.  'F'JTHEN 

RN-8 

ENDIF 

IF(CDL. EQ.  ' 0 ' >  THEN 

RN-13 

ENDIF 

IF(CDL.EO. 'H' >THEN 

RN-20 

ENDIF 

IF (CDL. EQ. 'J' JTHEN 

RN-32 

ENDIF 

IF(CDL.EO. 'K ' ) THEN 

RN-30 

ENDIF 

IF (CDL. EQ.  'L'JTHEN 

RN—80 

ENDIF 

IF( CDL.  EQ.  'fl '  J  THEN 

RN-123 

ENDIF 

IF(CDL.EO. 'N'JTHEN 
RN-20O  .  * 

ENDIF 

IF  (CDL.  EQ.  *1*  *  )THEN 

RN-3 IS 

ENDIF 

IFUCDL.  EQ.  'O'!  THEN 

RNtSOO 

ENDIF 

IF(CDL. EQ. 'R'lTHEN 

RN-800 

ENDIF 


IF(SP.EQ.  *D*. AND. CDL. EQ.  'A'JTHEN 

FS-0 

ENDIF 

IFISP.EQ.  'D '. AND. CDL. NE.  'A'JTHEN 

IF (CDL.  EQ.  'B'JTHEN 

FS-2 

ENDIF 

IF(CDL.EO.  'C'JTHEN 

FS-3 

ENDIF 

IF (CDL. EQ. 'D'JTHEN 

FS-3 

ENDIF 

IF(CDL.EO.  'E'JTHEN 

FS-8 

ENDIF 

IF(CDL. EO. 'F'JTHEN 

FS-13 

ENDIF 

IF (CDL. EO. 'C'JTHEN 

FS-20 

ENDIF 

IF(CDL. EQ. 'H'JTHEN 

FS-32 

ENDIF 

IF(CDL.EQ.  '  J '  1 THEN 

FS-30 

ENDIF 

IF(CDL.  EO.  'K'JTHEN 

FS-80 

ENDIF 

IF(CDL. EO.  'L '  J  THEN 


£ 


V 

V 


It 

,”1 


I 


'$ 


J 

'J 

U 


c 

c« 


K-2 

ELSE  IFtCDL.  EO.  'CM THEN 
Kv3 

ELSE  IFtCDL. EO.  *D')THEN 
K-4 

ELSE  IF(COL.EO.  *E')THEN 
K«3 

ELSE  IFtCDL.  EO.  *F*)THEN 
K-6 

ELSE  IF(CDL. EO.  'O') THEN 
K-7 

ELSE  IFfCDL.EO.  *H')THEN 

|^ag 

ELSE  IF(COL.EO.  'JMTHEN 
K«9 

ELSE  IF (COL. EO.  'K 'I THEN 
K»10 

ELSE  IFtCDL.  EO.  'L  *  ) THEN 
K=ll 

ELSE  IFtCDL. EO.  'M#> THEN 
K«12 

ELSE  IFtCDL.  EO.  'N'» THEN 
K»13 

ELSE  IFtCDL. EO. *P  ' ) THEN 
K«14 

ELSE  IFtCDL.EO. 'Q'lTHEN 
K-13 

ELSE  IFtCDL.  EO.  'R*)TH£N 
K«16 

ELSE  IFtCDL. EO.  *S*)THEN 

K-17 

END  IF 

IF  f AOL. EO. O. OIO  > THEN 
J»1 

ELSE  IFtAOL. EO.  0.  01 3) THEN 
J-2 

ELSE  IF t AOL. EO. . 023) THEN 
J-3 

ELSE  IFtAQL.EO.  0.040) THEN 
J-4 

ELSE  IFtAOL.  EO. O. 063>THEN 
J*>5 

ELSE  IFtAOL.  EO.  0.  10) THEN 
j-6 

ELSE  IFtAOL. EO. O. 13) THEN 
J-7 

ELSE  IFtAOL.  EO.  O.  23 ) THEN 
j-8 

EGSE  IFtAOL. EO. 0.4) THEN 
U«9 

ELSE  IFtAOL. EO. 0.63)THEN 
U*10 

ELSE  IFtAOL. EO. 1 >  THEN 
J-ll 

ELSE  IFtAQL.EO. 1.  5)THEN 
J-12 

ELSE  IF  t  AOL.  EO.  2. 5 ) THEN 
J-13 

ELSE  IFtAOL.  EO. 4) THEN 
J«14 

ELSE  IFtAOL.  EO. 6. 3) THEN 
J*15 

ELSE  IFtAQL.EO. 10) THEN 
U-16 

ELSE  IFtAQL.EO. 13)THEN 
J-17 

ELSE  IFtAQL.EO. 23) THEN 
U-18 

ELSE  IFtAOL.  EQ.40)THEN 
J-19 

ELSE  IFtAQL.EO. 631THEN 
J-20 

ELSE  IFtAOL. EO.  100) THEN 
J»21 
ELSE 
J«22 

END  IF 

RETURN 

END 

SUBROUTINE  VALUES 


K 

'i 
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FS-12S 
END  IF 

IF(CDL  EQ.'M' JTHEN 

FS-20L 

END  IF 

IF(CDL.  EQ.  'N  '  J THEN 

FS-313 

ENDIF 

IF(CDL.EQ.  'P' JTHEN 

FS-300 

ENDIF 

IF(CDL.EQ.  'Q' JTHEN 

FS-800 

ENDIF 

IF(CDL.EQ.  'R' JTHEN 

FS-12S0 

ENDIF 

IF(CDL.  EQ.  'S' JTHEN 

FS-2000 

ENDIF 

ENDIF 

IFUSP.EQ.  'D'J.AND. (CDL. Ed.  'B '. OR. CDL. EQ.  'C  *  >  JTHEN 
RFS-0 
ENDIF 

IF(SP.EQ.  'D ' . AND. CDL. NE.  'B AND. CDL. NE.  'C  '  JTHEN 
IF (CDL.  EQ.  'D' JTHEN 
RFS-2 
ENDIF 

IF(CDL. EQ.  'E  '  J  THEN 

RFS-3 

ENDIF 

IF(CDL.EQ.  'F' JTHEN 

RFS-3 

ENDIF 

IF (CDL. EQ.  'C' JTHEN 

RFS-8 

ENDIF 

IFiCDL.EQ. 'H' JTHEN 

RFS-13 

ENDIF  - 

IF (CDL. EQ.  'J' JTHEN 

RFS-20 

ENDIF 

IFCCDL.EQ.  'K' JTHEN 

RFS-32 

ENDIF 

IF(CDL. EQ.  'L' JTHEN 

RFS-30 

ENDIF 

IF (CDL.  EQ.  'M' JTHEN 

RFS-BO 

ENDIF 

IF (CDL. EQ. 'N' JTHEN 

RFS-123 

ENDIF 

IF(CDL.  EQ.  'P' JTHEN 

RFS-200 

ENDIF 

IF(CDL.  EQ.  'Q '  J  TJ(EN 

RFS-3 I 3 

ENDIF 

IF(CDL. EQ.  'R' JTHEN 

RFS-500 

ENDIF 

ENDIF 

RETURN 

END 

SUBROUTINE  INDEX 

**  THIS  SUBROUTINE  LOCATES  THE  CELL  IN  THE  *# 

**  TABLES  PRESENTED  IN  MIL-STD-103D  CONTAINC  ** 

**  ACCEPTANCE  AND  REJECTION  NUMBERS  FOR  ANY  ** 

**  PATICULAR  COMBINATION  OF  AQL  AND  CODE  LETTER.  *» 


CHARACTER  1*2. CDL.  SP 

INTEGER  FS. FA.  SA.  AR.  SR. AC. AT. AE. SC. ST, Z 
COMMON/SCL/L. SP. S.  I.  NQ.  T,  R,  CDL.  FS. D,  J.  K.  AQL. M.  FA.  SA. AR. 
*  SR,  AC.  AT.  AE.  SC,  ST.  Z 
IF(CDL. EQ. 'A 'JTHEN 
K-l 

ELSE  IF(CDL.EQ.  'B' JTHEN 
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C**  THIS  SUBROUTINE  FINDS  ACCEPTANCE  AND  REJECTION  ** 

C**  NUMBERS  FOR  THE  CELL  IDENTIFIED  BY  SUBROUTINE  INDEX.  •« 

C 

CHARACTER  1*2#  CDLi  SP 

INTEGER  FS. FA. SA. AR.  SR.  AC. AT.  AE.  SC.  ST.  Z.  NN.  NFS.  NFA.  NAE.  NSA 
INTEGER  NAR.  NAC. NAT.  TN.  TFS. TFA.  TAE.  TAR.  TAC.  TAT.  RN,  RFA.  RAE 
INTEGER  FRR.  FTR.  FNR.  SNR. STR, SRR.  RAR.  RAC.  RAT. RFS.  TSA 
COMMON/SCL/L.  SP.  S.  I.  NO.  T.  R.  CDL.  FS.  D.  J.  K.  AQL.  M.  FA.  SA. 

*  AR.  SR.  AC.  AT.  AE.  SC.  ST.  Z.  NTN.  NFS.  NFA.  NAE.  NSA.  NAR.  NAC.  NAT.  TN 

*  . TFS.  TFA. TAE. TSA, TAR, TAC. TAT. RN.  RFS.  RFA.  RAE.  RAR.  RAC.  RAT 

*  . FRR. FTR, FNR. SNR,  STR. SRR 
NFA- 100 

TFA- 100 
RFA- 100 
NAT-100 
TAT- 100 
RAT- 100 
TAC— 100 
RAC-100 
TAR- 100 
RAR -100 
TAE-lOO 
RAE-lOO 
TSA- 100 
NAC— 100 
NAR— 100 
NAE-lOO 
NSA— lOO 
ST-lOO 
SC-lOO 
SR-lOO 

IF1SP.E0. 'S') THEN 

IFIM.CE. 2. AND. M. LE. 161NFA-0 
IF (M. GE.  26 ) NFA— 21 
IF(M. EO. 23) NFA— 14 
IF (M. EQ. 24) NFA— 10 
IFIM.EQ.  23)NFA— 7 
IF  CM. EO.  22 ) NFA-3 
IF<M,.EQ.  21  )NFA— 3 
IFCM.EQ.  20) NFA— 2 
IF<M.  EO.  19)NFA-1 
/IFtM.EQ. 17. AND. J.  NE.  161NFA-0 
IFIM.EQ.  18.  AND. J. NE. 18)NFA-I 
>IF(K.  EQ.  16.  AND.  J.  LE.  21NFA-0 
lFlJ.EQ.  16. AND. K. LE.21NFA-1 
IF (M.  GE.  26) TFA- 18 
IFCM.EQ.  231TFA— 12 
~IF<M.  EQ.  24) TFA— 8 
_  IF<  M.  EQ.  23 ) TFA— 5 
V IF<M.  EO. 22)TFA— 3 
‘ IF<M.  EQ.  21 )TFA— 2 

IF(M. GE. 18. AND. M. LE. 20.  AND.  J.  NE.  161TFA-1 
IF (K.  EQ.  16.  AND. J. EQ. 2)TFA-2 
IFCM.LE. 17. AND. J. LE.  IS ) TFA— O 
IF<  J. EO.  16. AND.K. LE.  2) TFA- 1 
IF  <M. GE.  26)RAC— 10 
IF(M. EQ.  23)RAC— 7 
IFU1.EQ.  24)RAC-3 
IF CM. EQ.  23) RAC-3 
IF(M. EQ. 22.  AND. J. LE.  191RAC-2 
IF<M. EQ.  2DRAC-1 
IF (M. EQ. 20. AND. J. LE.  18)RAT-1 
IF(M. EO. 19.  AND. J. LE.  171RAT-0 
IF(M. GE. 2. AND. M.  LE.  16)RFA-0 
IF ( M. EQ. 17. AND. J.  LE.  1 5 ) RFA-0 
IF< J. LE. 2) RFA-0 
IF< J. EQ. 16.  AND.  K. LE.  31RAT-0 
IF (  J.  EQ.  17.  AND.  K.  EQ.  DRAT-0 
IF ( J. EQ. 2 1 .  AND. K. LE.  2 ) RFA- 5 
IFtJ.EQ.  20.  AND.  K.  EQ. l)RFA-3 
IF( J. EQ.  20.  AND.  K.  EQ.  2)RAT-3 
IF ( J.  EQ.  19.  AND.K.  EQ.  DRFA-2 
IF ( J. EQ. 19.  AND.  K. EQ.  2)RAT-2 
IF  ( J.  EQ.  18.  AND.K.  EQ.  DRFA-l 
IF ( J. EQ. 18. AND. K.  EQ.  2)RAT-1 

END  IF 

IFISP.EQ.  ' D ' >  THEN 

IF1M.LE. 16. AND. J. LE. 15) THEN 

NAT-0 

NSA-1 

TAT-0 

TSA-1 


RAT-O 
ST-0 
END  IF 

XFIK.Ed. 16. AND. J. LE. 2 ) THEN 
NAT-0 
NSA-1 
TAT-0 
TSA-1 
RAT-0 
ST-0 
ENDIF 

IFCK.EQ.  l.ANO. J. EG. IB 1 THEN 
NSA-3 
NAC-0 
TAT-0 
TSA-1 
END  IF 

IF<M.  CE.  26 > THEN 

NAE-11 

NSA-26 

END  IF 

IF  <11.  EQ.  23)  THEN 
NAR-7 
NS A— IS 
ENDIF 

IF(H.  EO.  24) THEN 
NAR-S 
NSA-1 2 
ENDIF 

IF<M. EO. 23) THEN 
NAR-3 
NSA-0 
ENDIF 

IFCM.EO.  22) THEN 
NAC-2 
NSA-6 
ENDIF 

IF<M.  EO.  21 )THEN 
NAC-1 
NSA-4  ' 

ENDIF 

XF<M. EO.  20 ) THEN 
.  NAC-0 
NSA-3 
V  ENDIF 

IF<  <N.  EQ.  IV.  AND. J.  CT. 2). OR. CM. EO. 18. AND.  J.  CT. 2) >THEN 
NAT-0 
NSA-1 
•1  ENDIF 

IFCJ.EQ. 16. AND.K.LE.31THEN 
v  NAT-0 
NSA-1 
ENDIF 

1F<M. CE. 26) THEN 
TAE-9 
TSA-23 
ENDIF 

IF<M. EQ. 29) THEN 
.  TAR-6 
TSA-1S 
ENDIF 

IF<M. EQ. 24) THEN 

TAR-3 

TSA-U 

ENDIF 

IF(M. EQ. 23) THEN 

TAC-2 

TSA-6 

ENDIF 

IF<M.  EQ.  22) THEN 
TAC-1 
TSA-4 
ENDIF 

IF<M.  EQ.  2DTHEN 


TAC-0 
TSA-3 
ENDIF 

IF  < M. CE. 18. AND.  M.  LE.  20.  AND. J.  CE.  3 ) THEN 
TAT-0 
TSA-1 
ENDIF 

IF<M.  EQ.  26) THEN 
RAE-3 
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I 


I 

I 


END  IF 
RETURN 
END 

SUBROUTINE  CODE 


SR- 12 
END  IF 

IF  CM. EO. 23) THEN 
RAE-3 
SR-B 
END  IF 

IF  CM. EO.  24) THEN 
RAE-2 
SC-6 
END  IF 

IF (M. EO. 23) THEN 
RAR-I 
SC-4 
END  IF 

IF  <M. EO. 22) THEN 
RAR-0 
SC-3 
END  IF 

IFCM.  EQ.2DTHEN 

RAR-0 

SR-1 

END  IF 

IF (M. EO. 20) THEN 

RAC-0 

SR-0 

END IF  : 

IFCM.  EO.  19) THEN 
RAT-O 
ST-O 

EIF C M.  EO.  IB.  AND.  J.  CE.  3.  AND.  J. LE. 1 3 ) THEN 
RAT-O 
ST-O 
ENDIF 


C****THIS*SUQROUT INE*IDENTIFIES  THE  CODE  LETTER  ASSOCIATED  ** 
C**  WITH  EITHER  THE  SPECIAL  OR  CENERAL  INSPECTION  LEVELS.  ** 
C**  FOR  ANY  GIVEN  LOT  OR  BATCH  SIZE.  _ _ _  ** 


CHARACTER  ’^*2.  CDL.  SP 
INTEGER  S1.S2.S3.S4 
COMMON/SCL/L.  SP.  S.  I.  NO.  T.  R,  CDL 
IFCI.E0~*S1 ')THEN 

IFCL. GE.  2. AND.  L.  LE.  30) CDL- 'A ' 

V  IFCL. GE.  31.  AND. L. LE. 3001CDL- 'B ' 
IFCL. CE. 301 . AND. L. LE. 33000 > CDL—  C 
IFCL. CE. 33001 >CDL-'D' 


ENDIF 

IFCI.EO. 


ENDIF 

IFCI.EO. 


'S2 ' )THEN 

IFCL. CE. 2. AND.  L.  LE.  23 > CDL-  A' 

IFCL. GE. 26. AND.L.LE.  130>CDL-'B' 
IFCL. GE. 131. AND.L.LE. 1200>CDL-'C' 
IFCL. GE. 1201.  AND.  L.  LE. 33000)CDL-'D‘ 
IFCL. GE. 33001) CDL- 'E' 


S3 ' ) THEN 

IFCL. CE. 2.  AND.L.LE.  13)CDL-'A' 

IFCL. CE. 26. AND. L.  LE.  30>CDL-'B ' 

I F  C L. CE. 3 1 .  AND.  L.  LE.  1 30 ) CDL- ' C ' 

IF  CL. CE. 1 31 . AND.  L.  LE.  300 >CDL- 'D ' 

IF  C L. CE. 501 . AND.  L.  LE.  3200 > CDL- 'E 
IFCL. CE. 3201.  AND.  L.  LE.  33000 ) CDL- 'F ' 
IFCL. CE. 33001. AND.L.LE. 300000) CDL- 'C' 
IFCL. CE. 300001 >CDL- 'H ' 

ENDIF 

IFCI.EO.  'S4  '  )THEN 

IFCL.CE.2. AND.L.LE.  13)CDL-'A' 

IF  CL. CE. 16. AND.L.LE.  23 >CDL-'B* 

IFCL. CE. 26. AND. L. LE. 90 )CDL- 'C * 

IFCL. CE. 91 .  AND.  L.  LE.  1 30 >CDL- 'D * 

IF  CL. CE. 131. AND.L.LE.  300)CDL-'E' 

IFCL. CE. 301.  AND.L.LE.  1200)CDL-'F' 
IFCL. CE. 1201.  AND.L.LE.  10000)CDL- 'C ' 
IFCL. CE. 10001. AND. L. LE. 3S000)CDL- *H' 
IFCL.CE. 33001. AND.L.LE. 300000>CDL-'J' 
IFCL. CT. 300000 ) CDL- 'K ' 

ENDIF 
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IFC I.EQ.  '  1  ' ) THEN 

IFCL.CE. 2.  AND.  L.LE. 13>CDL-'A* 

IF CL. CE. 16. AND. L. LE.  23) CDL” 'B ' 

IF < L. CE. 26. AND. L. LE. 90 ) CDL- ' C  ' 

IF (L. CE. 91 . AND. L. LE. 130)CDL-'D' 

IFCL.  CE.  131.  AND. L. LE. 280) CDL" 'E * 

IF(L.  CE.  281.  AND. L. LE. 300) CDL- *F  ' 

IF  (L.  CE.  301.  AND.  L.  LE.  1200  >CDL-  'C  * 

IF  < L.  CE.  1 20 1 .  AND.  L. LE.  3200 ) CDL- ' H  * 

IF (L.  CE.  3201.  AND.  L. LE.  10000)CDL-'U* 

IF  < L.  CE.  1 000 1 .  AND.  L. LE. 33000 ) CDL- ' K ' 

IF  (L.  CE.  33001 .  AND.  L.  LE.  1 30000  )CDL-  #L  * 

IFCL.CE.  130001.  AND. L.LE.  300000) CDL- *M* 

IF<L.  CE.  300001) COL- *N‘ 

END  IF 

IFCI.EO. *2') THEN 

IF (L. CE. 2. AND. L. LE.  8) CDL— 'A' 

IF(L. CE. 9. AND. L.LE. 13)CDL-*B* 

IF(L.  CE.  16. AND. L.LE.  23) CDL- 'C‘ 

IF ( L.  CE.  26. AND. L.  LE.  30 ) CDL- 'D ' 

IFCL. CE. 31. AND. L.LE.  90>CDL-'E' 

IFCL. CE.  91 , AND.  L.  LE.  1 30 ) CDL- 'F ' 

IF  <  L. CE. 1 31 . AND. L.  LE.  280)CDL- 'C ' 

IF ( L. CE.  281 . AND.  L.  LE.  300  >  CDL- *H ' 

IFCL.CE. 301. AND. L.LE.  1200)CDL-'J' 

IFCL. CE. 1201. AND. L.LE.  3200>CDL-'K* 

IF(L. CE. 3201. AND. L.LE.  10000)CDL- 'L ' 
IFIL.CE. 1000 1 . AND.  L.  LE.  33000 ) CDL- 'M ' 

I F ( L. CE. 3300 1 . AND.  L.  LE.  1 50000  >  CDL- *N  * 
IFCL.CE. 130001. AND. L.LE. 300000 > CDL- 'P * 

IFCL. CE. 300001 )CDL- 'G ' 

END  IF 

IFCI.EO.  '3  * ) THEN 

I F  C  L. CE. 2. AND. L. LE. 8 ) C  DL- *  B ' 

IFCL. CE. 9. AND. L.LE. 13>CDL-'C' 

IFCL.CE. 16. AND. L.  LE.  23 )CDL—  *D* 

IFCL. CE. 26. AND. L.LE. 30)CDL-'E  * 

IFCL.CE. 3 1. AND. L.LE. 90) CDL- #F * 

IFCL-.  CE.  91.  AND.  L.  LE.  130>CDL-'C' 

I F  C  L. CE. 1 3 1 . AND. L. LE.  280 ) CDL- #H ' 

IFCL. CE. 281. AND. L. LE.  300 > CDL- ' J ' 

IF C L.  CE. 30 1 . AND. L. LE. 1 200 ) CDL- ' K  * 

IFCL.  GE. 1201. AND. L.  LE.  3200) CDL- 'L ' 

I F  C  L. CE. 320 1 . AND. L. LE.  1 OOOO ) CDL- *  M ' 

I F  C  L. CE.  1 OOO 1 .  AND.  L.  LE.  3 3000 ) CDL- ' N  * 

IF  C  L. GE. 3300 1 . AND. L.  LE.  1 30000 ) CDL- 'P * 

IFCL. CE. 130001.  AND.  L.  LE.  300000 ) CDL- '0 ' 

IFCL.  CE.  300001 >CDL- 'R ' 

ENDIF 

RETURN 
ND 

SUBROUTINE  OCCNNNN.  IFS.  IRFS.  IRFA.  ISRR.  IFRR.NQ.  IRN.  ITFA. 
ITN.  IFTR.  ISTR.  NFA.  NTN.  IFNR.  I SNR.  IZ. J, AQL > 


THIS  SUBROUTINE  EVALUATES  THE  SCHEME  OPERATING 
CHARACTERISTICS.  FOR  THE  REQUIRED  SAMPLING  PLANS. 


REAL*8  AG C 20.  20). BO C 20. 20) 

DIMENSION  NC301.PRC30).  VC  100.  30).  PC30).  TPAC30.  4).  ASNC30.  4) 
DIMENSION  A03C30.  4).  ATIC30.  4).  NNC4).  TMC20.  20).  ZP C 20.  20) 
DIMENSION  ASNC  C 30 ) ,  AOQC  C 30 ) .  ATIC C 30 ) ,  AFI C 30 ) .  PAC 30 ) 

INTEGER  AC  30). RC 30).  HI.  AAC  2.  4).  RRC2.  4) 

IF C  J.  EG.  2)  CO  TO  197 

IFCJ.EQ. 1)  GO  TO  172 

NNC1 >-NQ 

NN  C  2 ) -NQ 

NNC3)-IRN 

NNC4 )— IRN 

AAC 1. 4)-IFRR—l 

RRC1.3)-IRFA+1 

CO  TO  272 

NNC 1 )-IFS 

NNC2)-IFS 

NNC3)-IRFS 

NNC  4 )  —  IRFS 

AAC 1 < 4 )-IRFA 

AA ( 2. 4  >  — ISRR— 1 

RRC 1. 3 ) — IFRR 

AAC 1 . 1)-ITFA 

AAC 2. 1 >— ITN 

RRC 1. 1 )-IFTR 


v.vw, 


uLMjuuu  uuuuuu  uuuuuu 


RRC2. D-ISTR 
AAC1.21-NFA 
AA(2.  2>-NTN 
RRU.  2>-IFNR 
RR  <  2. 2J-ISNR 
AA< i.  3>*IRFA 
AA<2.  3)-I  Z 
RR ( 2. 3>-IZ+l 
RR< 1> 4>— IFRR 
RRC2.  4>— ISRR 

PRINT  *. 'DO  YOU  WANT  A  TABLE  OR  A  GRAPH  FORMAT  ?' 

PRINT  *. '  FOR  GRAPH. ...ENTER:  1  ' 

PRINT  *.  *  FOR  TABLE. ...ENTER:  2  ' 

READ  •. XTC 

IF <  XTC. EQ. 1 )  GO  TO  733 

PRINT  *.  'SPECIFY  THE  NUMBER  OF  FRACTION  DEFECTIVE  VALUES.. 
READ  *.  J1 

PRINT  *.  'ENTER  THE  FRACTION  DEFECTIVE  VALUE(S).  ' 

PRINT  *.  '(PUT  A  COMMA  BETWEEN  VALUES.  > . ' 

READ  *. <P< I >. I-X. J1 > 

WRITE  <3. 27) 

CO  TO  447 
Jl-21 

DO  33  I-l.Jl 
P(I)-(I-l>/20.0 
P(I)-P(X>/10.0 
DO  103  L-l. J1 
DO  tO  LM-1.4 
DO  17  I-l.J 
AC  I ) «AA( I •  LM ) 

RCI)-RR(I.LM) 

NCI J-NNCLM) 

CONTINUE 

NNN-NNNN 

H-INT(NNN*PCL>> 

ZNK-NNN*PCL> 

ZN-M 

IFC  ZNK-ZN. LE. 0.  3)  CO  TO  148 

M-M+l 

CONTINUE 

SUM-O 

Il-RCl) 


THE'  NEXT'  STATEMENTS'  CALCULATE' THE'  PROBABILITY  THAT' 
THERE  ARE  <I-t>  DEFECTIVES  IN  THE  FIRST  SAMPLE. 


DO  20  I-l. It 

V(I.  1 )— PPCNC 1 ) i  P (L) •  I-l) 
SUM-SUM+VI I.  1 } 

VCRC  1)  +  I.  D-X.-SUM 
IFCACD.LE.  0)  GO  TO  2X 
Il-ACX) 


THE  FOLLOWING  STATMENT8  CALCULATE  THE  PROB.  THAT  THERE 
ARE  LESS  THAN  OR  EQUAL  TO  A(X>  DEFECTIVE  ITEMS. 


DO  30  I-l. IX 


THE  NEXT  SECTION  CALCULATES  THE  INTERMEDIATE  PROBABILITIES 
OF  CONTINUED  SAMPLINO  FOR  THE  DOUBLE  SAMPLING  PLANS. 


V(A(  1  )♦!.  1)>V(A<1HXi  1)+VCI.  X) 
IFC J. EO. X)  CO  TO  41 
LOW  *1 
1-2 

NNN-NNN-NC I-l > 

NR-R ( I > 

LOW-MAX 0 ( LOW.  A  C I - 1 >  +2 ) 

HI -R (1—1) 

DO  60  11-1.  NR 
VCIl, I>-0 

IFCI1.LT. LOW)  GO  TO  60 
IFCACD.EQ.  C  — 1  >  >  GO  TO  61 
IFCACIH-l.NE.  ID  GO  TO  61 
MMM-MINOC 1 1.  HI ) 

DO  70  12-LOW  .  MMM 


I4-I1-I2+1 
DO  SO  II-l.NR 

PR!  II)-PP(N(I>.P!L).  II-l.NNN.  M-12-M) 

SO  CONTINUE 

DO  80  13-1. 14 

80  VIII.  I)«V!ll,  I)+V!I2.  I— 1 )  *PR  ( 13) 

70  CONTINUE 
CO  TO  60 

61  MMM-MINO(HI,  Ill 

DO  90  12-LOW  <  MMM 

I3-I 1-I2+1 
DO  743  IJ-l.NR 

743  PRCIJ)«PP!N(I)«P(L>.  IJ-1) 

90  VIII.  I)-V(I1.  I)+V!I2.  I  — 1 ) *PR ( 13 ) 

60  CONTINUE 

V(RCI)-»-l.  I)-0 
DO  110  Il-LOU  .HI 
SUM-0 

I3-R ( I )— Il+l 
DO  120  12-1. 13 

120  SUM-SUM+PR(I2> 

IlO  V(R( 1 1  +  1.  1 )-V!R< I )+l*  I )  +  !  1.  -SUM)*V( II.  I-l> 

41  CONTINUE 
C 

C:  :::::::::::::::::::  :  ::::::::::::::::::::::::::::::::::  : 

C  THE  NEXT  SECTION  EVALUATES  THE  SCHEME  FOR  THE  SINGLE 
C  SAMPLINC  PLAN  (OR  THE  FIRST  STACE  IN  DOUBLE  SAMPLING). 

C:  ::::::::::::::::::::::::::::::::::::::::::::::::::::::: : 

c 

121  ASN(L. LM)-0 
TPA(L. LM>— O 
SS-O 

DO  130  I-l.J 
SS— SS+N  < I > 

IF(A(I).EQ. <— 1>>  GO  TO  131 

ASN  <  L.  LM  )  — ASN  (  L.  LM)«-!V!AC  I  »♦!.  I  >*V(R(  I  )♦!.  I) >*SS 
TPACL.  LM)-TPA!L.  LM)«-V!A(I)»1,  I) 

GO  TO  130 

131  ASN(L. LM) — ASN(L»  LM)4V(R ( I  )  +  l .  I)*SS 

130  CONTINUE 
XXX-O 
SSS-O 
WV-0 

DO  133  I-l.J 
SSS-SSS4-/K  I ) 

XXX-XXX+SSS*VCA( I )♦!• I) 

WV-WV+  ( NNNN-SSS  )  *V  <  A  <  I  )  ♦  1 ,  I) 

133  CONTINUE 

ATI  CL. LM)— XXX-MNNNN*! 1 . 0—TPA(L.  LM) ) > 

AOO(L<  LM)  — (WV*P(L)  l/NNNN 
140  CONTINUE 
lO  CONTINUE 

19  FORMAT! IX.  F3. 3.  2X.  F9.  4.  FlO. 2.  4X.  F6. 4.  IX.  FIO. 2) 

27  FORMAT!///'  SCHEME  OPERATING  CHARACTERISTICS 

:P(A>: : : : : : :  ASN:  : : : : : AOQ: :::::::: :  AF I ' > 

DO  1  1-1.20 
DO  1  JJ-1.20 
1  TM( I. JJ)-0 


C  THIS  SECTION  FIRST  CALCULATES  THE  PROB.  OF  BEING 
C  IN  NORMAL.  TIGHTENED  OR  REDUCED  INSPECTION  LEVELS  AND 
C  THEN  COMPUTES  THE  PROB.  OF  PASSING  THE  LIMIT  NUMBER 
C  CRITERIA. 


T-TPACL. 1 ) 

S-TPACL. 2) 

RRl-TPACL.  3) 
Rl-TPACL.  4) 
NZ“10*ASN(L.  2) 

CALL  LNMCAQL.  NZ.LNC) 
Q-0 (NZ.P(L).  LNC ) 


THIS  SECTION  READS  IN  THE  PROBABILITY  TRANSITION  MATRIX. 

TM! 1.1 )-l.  -T 
TM( 1. 2)-T 
TM ( 2. 1 )-l. -T 
TM(2.  3)-T 
TM(3.  1 )-l. -T 


HPB 


U<JO<JU<JU 


TM<3. 4)-T 
TM<4.  1 )-l. -T 
TM<4.  3)-T 
THIS.  D-l.-T 
TH<3.  6>«T 
TM<6.  7>-l.  -S 
Tt1(6.  1D-S 
TH<7.  D-l.-S 
TM<7.  8>»S 
TM<8. 1>»1.-S 
TM<8.  9>-S 
TM<9.  1 >-l. -S 
TM(9.  10>-S 
TM( 10.  D-l.-S 
TM( 10#  14  >— S 
TM (11# 7)— 1. — S 
TM< 11# 12) — S 
TM< 12#  7)  — 1.  -S 
TM< 12. 13)— S 
TM< 13.  7>-l.-S 
TM(  13.  14)-S 
TM< 14.  7)  — 1.  -S 
TM(14.  13)-S 
TM( 13.  7)-l.  -S 
TM< 13.  16)— S 
TM(  16.  7)-l.  -S 
TM< 16. 17)-S 
TM<17,  7)-l.-S 
TMD7.  18)— S 
TM<  18.  7>-l.  -S 
TM(1B. 19)-S 
TM< 19, 7)— 1.  — S 
TM< 19. 19>— ( 1. -0)*S 
TM< 19. 20 ) -Q*S 
TCI <20.  6)-l.  -RR1 
TM( 20. 20)“RR1 
NR— 20 
NC-20 

DO  100  I-l.NR 
DO  lOO  UJ-l.NC 
IF(JJ-I)111. 13. Ill 
ZP(JJ.  D-TMU.  JJ)— 1.0 
GO  TO  lOO 
ZP<UJ.  M-THU,  JJ) 
CONTINUE. 

DO  113  r^2.NR 
DO  113  JJ-l.NC 
KL-I-1 

AQ< I.  JJ)-ZP<KL.  JJ) 
CONTINUE 
DO  18  »JJ— 1 .  NC 
AQ<  1 ,  JJ)-1.0 
CONTINUE 

CALL  INVERT1NR.  AO.  BQ) 

PA(L)-0 

ASNC  <L  >  — O 

AOQC ( L ) — O 

ATIC (L)— O 

DO  161  I-l.NR 

IF<  I.CT.'3)C0  TO  200 

ZZ-T 

Z1-ASN(L,  1) 

Z2— A00(L. 1 ) 

Z3-ATI <L.  1) 

GO  TO  201 

IF < I.GT. 19)  GO  TO  210 
ZZ-S 

Z1-ASN1L. 2) 

Z2-A00(L. 2) 

Z3-ATKL.  2) 

CO  TO  201 
ZZ-R1 

Zl-ASNCL. 4) 

Z2-AQ0(L«  4  > 

Z3-ATI <L.  4) 


THE  FOLLOWING  STATEMENTS  MULTIPLY  THE  PROBABILITIES 
OF  BEING  IN  THE  VARIOUS  STATES  BY  THE  CHARACTERISTICS 
OF  EACH  STATE. 


PA(L)-PA(L)«-BO(  I.  1  )*ZZ 


59 


,TB 

r? 


I 

$ 

S 

: s 


S 

J 


I 

« 

6 


ASNC ( L ) —ASNC ( L >  +B0 ( I .  1>*ZI 
AOQC!L)-AOQC!L>«-BQ! I,  1)*Z2 
AT I C !  L ) “AT I C ! L ) +BQ ! I » 1)*Z3 
AFI!L)-ATIC!L)/NNNN 
161  CONTINUE 

IF1XTC.NE.  2)  CO  TO  109 

WRITE (3. 19)  P  <L)  •  PA(L  >  •  ASNC(L) • AOOC (L)  •  AFI <L) 
105  CONTINUE 

IF ( XTC. EQ. 2)  CO  TO  596 
WRITE13. 45) 

WRITEIS. 46) 

CAUL  PLOTOC IP.  PA) 

WRITE t 5. 47) 

CALL  PLOTOC !P.  ASNC) 

WRITE! 5. 48) 

CALL  PLOTOC IP.  AOQC) 

WRITE (3.  49) 

CALL  PLOTOC (P»  AFI ) 

45  FORMAT! 1H1.  SOX.  'OPERATING  CHARACTERISTIC  CURVE*) 

46  FORMAT! 53X.  'PROBABILITY  OF  ACCEPTANCE*) 

47  FORMAT! 1H1.  SOX.  'AVERAGE  SAMPLE  NUMBER') 

48  FORMAT! 1H1.  SOX.  'AVERAGE  OUTGOING  QUALITY') 

49  FORMAT!! HI.  SOX.  'AVERAGE  FRACTION  INSPECTED') 

956  CONTINUE 

RETURN 

END 

FUNCTION  PP!N.  P.  K) 

0*1.  -P 
PP*Q**N 

IF !K.EO.O> RETURN 
DO  20  1*1 »  K 

20  PP*PP*P*!N-I+1 )/!0*I > 

RETURN 

END 

FUNCTION  BIN. P. K> 

PN*P*N 
0-1.  -P 

IF!PN. CT.  3.  )  GO  TO  1 
B-EXP1-PN) 

Z-EXPI-PN) 

IF IK. EQ. O)  RETURN 
DO  lO  1-1. K 

z— z*pn/j 

lO  B-B+Z 

RETURN  . 

1  B— Q**N 

Z— Q**N 

IF IK.EQ.O) RETURN 
DO  20  I-— 1 «  K 
Z*Z*!N-I+1 )*P/!Q*I > 

20  B-B+Z  < 

RETURN 

END 

SUBROUTINE  INVERTINR.  AQ. BQ) 


C** 

C** 


THIS  SUBROUTINE  INVERTS  THE  MATRIX  BY  GAUSSIAN 
ELIMINATION. 


** 

** 


lO 

is 


20 


29 

30 
33 


C 

c* 


REAL*8  AQ! 20.  20).  BQI20.  20).  ZTAMP. ATEMP 

DO  IS  1*1. NR 

DO  10  J-l.NR 

BO!I. JI-0.0 

BQ!  I >  D-l.O 

DO  3S  I-l.NR 

ZTAMP—AQ! I .  I) 

DO  20  J-l.NR 

AQ! I, J>— AQ! I.  J) /ZTAMP 

BQ! I. J)-BQ! I.  J) /ZTAMP 

DO  30  I I-l.NR 

IF! I. EO. II)  GO  TO  30 

ATEMP-AQ III.  I) 

DO  29  J-l.NR 

AO! II. Jl-AQIII.  J)-AQ!I.  J)*ATEMP 

BQ!II.J)-BQ!II.J) — BQ ! I .  J  >  * ATEMP 

CONTINUE 

CONTINUE 

RETURN 

END 

SUBROUTINE  PLOTOC !P.T) 
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'jrv  /  / 


Sl'./.WA. 


THIS  SUBROUTINE  PLOTS  OC.ASN. AOQ  AND  AFI  CURVES  FOR  * 
THE  FULL  RANGE  OF  INCOMING  FRACTION  DEFECTIVE.  « 

»*****#************##**#********************************< 

DIMENSION  U(  ID.  TOO).  POO).  LINE!  101) 

INTEGER  ASTERK.  BLANK.  PLUS 

DATA  ASTERK. BLANK.  PLUS/'*'.  '  '.'♦'/ 

DO  1  K-l.  101 
LINE (K) -BLANK 
XMAX-T  C 1 ) 

XMIN-T ( 1 ) 

DO  2  1-1.21 

IFCTCI). LT. XMIN)  XMIN-T ( I ) 

IFCTCI). CT.XMAX)  XMAX-TU) 

CONTINUE 

IFIXMAX.LT. 1 >  XMIN— O 
RANCE-XMAX— XMIN 
RC— RANGE/10 
U( 1 )— XMIN 
DO  3  1-2. 11 
UCI)-UCI-1)+RG 

WRITEO.  130)(U(1).  1-1.  11)  . 

DO  4  K-l. 101. 10 
L I NE  <  K  > -PLUS 

WRITEO.  100)  (LINE(K).  K-l.  101  > 

DO  6  K-l, 21 
DO  3  1-1, 101 
LINE! I >— BLANK 

KPRINT— lOO* <T(K) -XMIN) /RANGE* 1. 3 
LINE! 1 ) —PLUS 
L I NE  <  KPR I NT ) -ASTERK 

WRITEO.  1 40)P <K > .  (LINED).  1-1.  101) 

CONTINUE 

FORMAT (4X.  'PERCENT  DEFECTIVE '.  101A1 ) 

FORMAT ( 12X.  1 1  <2X. E8.  3) ) 

FORMAT ( 17X. F4. 3.  101A1) 

RETURN 

END 

SUBROUTINE  LNMCAQL.  NZ.  LNC) 

THIS  SUBROUTINE  FINDS  THE  LIMIT  NUMBER  ** 

FOR  REDUCED  INSPECTION.  ** 

V 

LNC-222 

IF(NZ.LE. 29. AND. AQL.LE.  13)  LNC-0 
IF<NZ.GE.30. AND.NZ.LE.  49. AND. AOL.LE.  10)  LNC-0 
IF(NZ. CE.  30. AND.  NZ. LE.  79.  AND.  AQL. LE. 6. 3)  LNC-O 
IFCNZ. GE. 80. AND.NZ.LE.  129.  AND.  AQL. LE. 4. O )  LNC-0 
IFCNZ. GE. 130. AND.NZ.LE.  199.  AND. AQL.  LE. 2. 3)  LNC-O 
IF < NZ. GE. 200. AND.  NZ.  LE.  319.  AND.  AQL.  LE. 1.  3)  LNC-O 
IFCNZ. GE.320.  AND.NZ.LE.  499. AND. AQL.LE. 1.0)  LNC-O 
IFCNZ. GE. 300. AND.  NZ.  LE.  799.  AND. AQL.  LE. 0. 65)  LNC-0 
IFCNZ. GE. 800. AND.  NZ.  LE.  1249.  AND.  AQL.  LE.  0.  40)  LNC-0 
IFCNZ.  GE.  1230.  AND.NZ.LE.  1999.  AND.  AQL. LE. 0. 23)LNC«0 
IFCNZ. GE. 2000.  AND.  NZ.  LE.  3149.  AND.  AQL.  LE.  O. 13)LNC-0 
IFCNZ.  GE.  3130.  AND.  NZ.  LE.  4999.  AND.  AQL. LE. O. 101LNC-0 
IFCNZ.  GE.-SOOO.  AND.  NZ.  LE.  7999.  AND.  AQL.  LE.  0.  0631LNC-0 
IFCNZ. GE.  8000.  AND.  NZ.  LE.  12499.  AND.  AQL. LE. O. 0401LNC-0 
IFCNZ.  GE. 12300.  AND.  NZ. LE.  19999.  AND.  AQL.  LE.  O.  025)LNC-0 
IFCNZ. GE. 20000.  AND.  NZ.  LE.  31499.  AND.  AQL.  LE. 0. 01 3) LNC-0 
IFCNZ. GE. 31300.  AND.  NZ.  LE.  49999.  AND.  AQL.  EQ.  0. 010)LNC-0 
IFCLNC.EQ.  O)  CO  TO  lO 
IFCNZ. LE. 29) THEN 

IFCAQL. EQ. 25>LNC-2 
IF C AQL.  EQ.  401LNC-4 
IFCAQL. EQ.  63 )LNC-8 
IFCAQL. EQ.  100)LNC— 14 

END  IF 

IF  C NZ. GE. 30. AND.  NZ .  LE.  49 ) THEN 
IFCAQL. EQ. 1 5 )LNC— 1 
IFCAQL.  EQ.  25)LNC-3 
IFCAQL. EQ. 40>LNC— 7 
IFCAQL. EQ.  63 1LNC-13 
IFCAQL. EQ.  100 )LNC-22 

END  IF 

IFCNZ. GE. 30. AND.  NZ.  LE.  79  > THEN 
IFCAQL. EQ.  101LNC-2 
IFCAQL. EQ.  13) LNC— 3 
IFCAQL. EQ.  23>LNC-7 
IFCAQL. EQ.  40 >LNC  —  1 4 
IFCAQL. EQ.  651LNC-23 


IF<AQL. EO. 100)LNC-40 

END  IF 

IFCNZ. CE. 80. AND.NZ.LE.  129) THEN 
I F  <  AOL. EO.  6.  9 ) LNC-2 
IFCAQL.EQ.  10>LNC-4 
IF (AOL. EQ. 19>LNC-7 
IF  (AOL.  EQ.29)LNC-14 
IF (AOL. EQ.40ILNC-24 
IF  ( AOL.  EO.  69  >LN042 
IF ( AOL.  EO.  100)LN068 

END  IF 

IFCNZ.CE. 130. AND.NZ.LE.  199) THEN 
IF  (AOL.  EO.  4.  0>LNO2 
IFCAQL.EQ. 6.  3) LNC-4 
IF (AOL. EO. 10ILNC-7 
IFCAQL.EQ.  15)LNC-13 
IF(AQL.EQ. 29ILNC-29 
I F ( AQL. EQ.  40  > LNC-4 2 
IF  (AOL.  EO.  691LN072 
I F (  AOL.  EO.  100ILN0119 

END  IF 

!F(NZ. CE. 200.  AND.  NZ.  LE.  319) THEN 
IF  (AQL.  EQ.  2.  9  )LNC*>2 
IF ( AQL. EQ.  4.0) LNC-4 
IFCAQL.EQ.  6.  9>LNC-B 
IF (AOL. EO. 10ILNC-14 
IF ( AQL. EQ: 13ILNC-22 
IF ( AQL. EO. 29)LNC-40 
IFCAQL.EQ.  40ILNC-68 
IF  (AOL,  EQ.  691LN01 19 
IF  (AQL.  EO.  100)LN0181 

ENDIF 

IF (NZ. CEi 320.  AND.  NZ.  LE.  499 > THEN 
IFCAQL. EQ.  1.9) LNC  * 1 
IF (AOL.  EO.  2.  5)LNC=4 
IF (AOL.  EO.  4>LNC“8 
I F  (  AQL.  EQ.  6.  9  )  LNC- 1 4 
IF  (AOL.  EQ.  10JLNO24 
IFCAQL.EQ.  13ILNC-39 
I F ( AQL.  EQ.  29 1 LNC-68 
IF (AQL. EQ. 631LNC-U.3 
IFCAQL.EQ.  100)LNC-189 

ENDIF 

IFCNZ. CE; 900. AND.  NZ.  LE.  799  > THEN 
IFCAQL. EQ. 1.0)LNC>2 
IFCAQL.EQ.  1.  91LNC-3 
IFCAQL. EQ. 2. 9JLNC-7 
IFCAQL.EQ.  4.0ILNC-14 
IF.CAQL.  EQ.  6.  9)LNC-29 
IFCAQL.EQ. 10ILNC-40 
IFCAQL.EQ.  19)LNC-63 
IFCAQL. EQ. 29) LNC* 1 10 
IFCAQL.EQ.  40ILNC-181 

ENDIF 

IFCNZ. CE. 800. AND. NZ . LE.  1 24  9  >  THEN 
IF  C AQL. EQ. 0. 69 ) LNC-2 
IFCAQL.EQ. l.OILNC-4 
IFCAQL.EQ. 1. 9)LNC-7 
IFCAQL.EQ.  2. 9>LNC«14 
IF'C  AQL.  EQ.  4.  O )  LNC-24 
IFCAQL. EQ.6. 9)LNC-42 
IFCAQL.EQ. 101LNC-68 
IFCAQL.EQ. 19)LNC-109 
IFCAQL. E0.291LNC-1 81 

ENDIF 

IFCNZ. CE.  1290. AND.NZ.LE.  1999) THEN 
IFCAQL.EQ.  0.  40ILNC-2 
IFCAQL. EQ.  0. 69)LNC-4 
IFCAQL.EQ. 1 . 0 ) LNC *7 
IFCAQL.EQ. 1. 9ILNC-13 
IFCAQL. EQ. 2. 91LNC-24 
I F ( AQL. EQ. 4. O  > LNC-40 
IFCAQL.EQ.  6.  9) LNC *69 
IFCAQL.EQ. 10)LNC-110 
IFCAQL.EQ. 19ILNC-169 
ENDIF 

IFCNZ. CE. 2000. AND. NZ. LE. 3149 1THEN 
IFCAQL.  EQ. 0.  29)LNC«2 
IFCAQL. EQ.  0. 40) LNC-4 
IFCAQL.  EQ. O. 69)LNC»8 
IFCAQL.EQ.  1 . 0  > LNC  —  1 4 
IFCAQL.EQ. 1. 9ILNC-22 
IFCAQL. EQ. 2. 9)LNC=40 


IF( AQL. EQ. 4. 0)LNC«68 
IF( AQL. EO. 6. 5) LNC- 103 
IF( AQL. EQ.  101LNC-181 

END  IF 

IF (NZ. CE. 3130. AND. NZ. LE. 4999 1THEN 
IF (AQL. £0.  0. 13JLNC-1 
IF ( AOL. EQ.  0.  23  > LNC-4 
IF (AQL. EQ. O. 40>LNC-8 
IF (AQL. EQ.  0. 631LNC-14 
I F  ( AQL.  EQ.  1 .  O  >  LNC-24 
IF  (AQL .  EQ.  1.  31LNC-38 
IF (AQL. EQ. 2. 31LNC-67 
IF ( AQL. EQ. 4. 0)LNC-1 1 1 
IF( AQL. EQ.  6. 5)LNC-186 

END  IF 

IF(NZ. GE. 3000. AND.  NZ.  LE.  7999) THEN 
IF  (AQL.  EQ.  0.  10)LNO2 
IF  (AQL.  EQ.  O.  13>LN03 
IF (AQL. EQ. 0. 23>LNC»7 
IF (AQL. EQ. 0. 40) LNC- 14 
IF( AQL. EQ. 0. 631LNC-25 
IF(AQL. EQ.  1.0)LNC-40 
IF ( AQL. EQ.  1.  5 ) LNC  =63 
IF( AQL.  EQ.  2.  3)LNC-1 10 
IF (AQL. EQ.  4.  O) LNC- 181 

END  IF 

IF ( NZ. GE. 8000. AND. NZ.  LE.  12499) THEN 
IF( AQL. EQ. O. 063 )LNC— 2 
IF ( AQL. EQ. O. 10 ) LNC— 4 
IF( AQL. EQ. 0. 151LNC-7 
IF(AQL.EQ.O.  25>LNC*-14 
IF ( AQL. EQ.  O.  40 ) LNC-24 
IF (AQL. EQ.  O.  63 ) LNC-42 
IF ( AQL. EQ. 1.0)LNC«68 
IF(AQL. EQ.  1.5) LNC* 105 
IF ( AQL. EQ. 2. 31LNC-181 

END  IF 

IF(NZ. GE. 12300.  AND.  NZ.LE.  199991THEN 
IF (AQL.  EQ.  0.  040 )LNC=2 
IF(AQL.  EQ.  O.  063) LNC-4 
IF(AQL. EQ.  O.  10)LNC«7 
IF (AQL. EO. O. 13) LNC* 1 3 
IF( AQL. EQ. 0.  2S)LNC«24 
IF ( AQL.  EQ. O. 40 ) LNC-40 
.  IF (AQL. EQ. O. 63 ) LNC— 69 
•  IF (AQL. EQ. 1. 0)LNC-1 10 
IF< AQL.  EQ.  l.3)LNC«169 

ENDIF 

IF (NZ. GE. 20000. AND.  NZ.  LE.  314991THEN 
IF (AQL. EQ. 0.  023 )LNC-2 
IF( AQL. EQ. 0. 040)LNC=4 
I F ( AQL. EQ. O. 06  3  >  LNC-8 
IF( AQL. EQ. 0. 10 ) LNC— 1 4 
IF(AQL. EQ.O.  13)LNC=22 
IF (AQL. EQ.  O.  25 )LNC*40 
IF( AQL. EQ. O.  40>LNC*6B 
IF  (AQL.  EQ.  0.  63  )LNC*>1  1  5 
IF( AQL. EQ. 1. 0)LNC“181 

ENDIF 

IF(NZ. GE. 31500. AND.  NZ. LE.  49999) THEN 
IF (AQL. EQ. 0. 01 3) LNC- 1 
IF ( AQL. EQ. 0. 025 > LNC-4 
IF ( AQL. EQ. O. 040 ) LNC-8 
IF(AQL. EQ. O. 063  >LNC— 14 
IF(AQL. EQ,  0. 10 > LNC— 24 
IF( AQL. EQ. 0. 1 3 ) LNC-38 
IF (AQL. EQ. 0. 25 )LNC=67 
IF (AQL. EQ. 0. 40) LNC- 1 1 1 
IF ( AQL. EQ. 0. 65) LNC- 186 

ENDIF 

IF ( NZ. GE.  50000) THEN 

IF (AQL. EQ.O. 010) LNC— 2 
IF (AQL. EQ. 0. 015 )LNC-3 
IF( AOL. EQ.O. 025 >LNC=7 
IF(AQL. EQ. 0. 040) LNC- 14 
IF (AQL. EQ. 0. 063 ) LNC— 23 
IF ( AQL. EQ. 0. 10 ) LNC— 40 
IF ( AQL. EQ. 0. 15) LNC— 63 
IF (AQL. EQ. 0. 25 ) LNC- 1 1 0 
IF (AQL. EQ. 0. 40  >  LNC— 181 
IF( AQL. EQ. O. 63 ) LNC— 30 1 

ENDIF 

CONTINUE 

RETURN 

END 


1 


